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boundary-layer  characteristics;  streamwlse  turbulence 
Intensities;  and  fluctuating  surface  pressures  have  been 
spatially  mapped  In  the  region  of  each  edge.  Correspond- 
ence between  surface  pressures  and  turbulent  fields  In 
both  the  Immediate  vicinity  and  the  far  wake  of  each  edge 
were  determined  by  correlation  analyses.  Results  show  a 
mechanism  of  pressure  production  which  Is  common  to  both 
random  and  periodic  pressures.  They  also  Indicate  that 
the  mechanism  differs  from  that  which  produce^  attached 
boundaiy-layer — wall-pressure  fluctuations. <^he  level 
and  chordwlse  variations  of  periodic  pressure^v^e  shown 
to  be  determined,  through  a nearly  unique  fuactlohal 
form,  by  the  strength  and  formation  length  of  vortlc^a^ 
Furthermore,  the  vortex  shedding  frequency  Is  shown  to  ae 
expressible  In  alternative  dimensionless  forms  based  on 
elthe:  near-  or  far-wake  variables.  In  the  cases  of 
random  wake  generation,  the  surface  pressures  are  generated 
by  both  a large-scale  eddy  structure  and  a smaller  scale 
convected  field.  The  large  eddies  are  formed  In  the 
separation  process  at  the  edge;  but  they  do  not  form  the 
coherent  wavelike  wake  structure  which  Is  shown  to  be 
necessary  for  the  generation  of  periodic  pressures.  Al- 
though the  broadband  Intensities  of  the  turbulence  In 
the  near  wakes  of  both  types  of  edge  have  been  found  to 
be  comparable,  a coherent  wake  structure  Is  seen  to  exist 
behind  only  those  edges  sustaining  periodic  pressures. 
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ADMINISTRATIVE  INFORMATION 


The  work  reported  herein  was  funded  by  the  Naval  Sea  Systems  Command 
under  Subproject  SF43-452-702 , Task  18184,  Work  Unit  1-1942-080  and  by 
the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center  In- 
House  Research  Program  under  Subproject  ZR  Oil  08  01,  Work  Unit  1-1942-069. 

1.  INTRODUCTION 

Influences  of  trailing-edge  shapes  on  vortex-induced  vibrations  of 


turbine  blades  were  the  subjects  of  an  early  investigation  by  Heskestad 

and  Olberts  (1960).^  In  that  series  of  measurements  it  was  found  that  the 

most  severe  hydroelastic  vibrations  of  a flat  plate  were  Induced  by 

blunted  trailing  edges  and  edges  symmetrically  beveled  or  tapered  with 

large  Included  angles  at  the  apices.  The  lowest  levels  of  vibration 

occurred  for  edges  which  were  unsymmetrically  beveled  with  relatively 

small  included  angles  and  a large  curvature  at  the  juncture  of  the 

beveled  edge  with  the  main  surface  of  the  strut.  Toebs  and  Eagleson 
2 

(1961)  not  only  characterized  some  effects  of  trailing-edge  geometry  on 
vibrations  but  also  experimentally  demonstrated  some  aspects  of  self- 
excitation  in  the  hydroelastic  behavior  of  flat  plates.  Eagleson 

3 

et  al.  (1962)  related  those  results  to  the  development  of  periodic 
disturbances  in  the  wake,  finding  that  the  strengths  of  shed  vortices 
could  be  reduced  by  increasing  the  taper  of  the  edge.  They  also  found 
that  if  the  plate  were  allowed  to  vibrate,  the  amplitude  of  motion  in- 
duced by  trailing-edge  flow  would  decrease  commensurately  with  the  ob- 
served reduction  in  the  strength  of  shed  vortices. 


■^^Heskestad,  G.  and  D.R.  Olberts,  "Influence  of  Trailing  Edge  Geometry 
on  Hydraulic-Turblne-Blade  Vibration  Resulting  from  Vortex  Excitation," 
Transactions  of  the  American  Society  of  Mechanical  Engineers,  Journal  of 
Engineering  For  Power,  Vol.  82,  pp.  103-110  (Apr  1960).  A complete  listing 
of  references  is  given  on  page(s)  152. 

2 

Toebs,  G.H.  and  P.S.  Eagleson,  "Hydroelastic  Vibrations  of  Flat  Plates 
Related  to  Trailing  Edge  Geometry,"  Transactions  of  the  American  Society  of 
Mechanical  Engineers,  Journal  of  Basic  Engineering,  pp.  1-8  (1961). 

3 

Eagleson,  P.S.  et  al.,  "Turbulence  In  the  Early  Wake  of  a Fixed  Flat 
Plate,"  Massachusetts  Institute  of  Technology  Hydrodynamics  Laboratory 
Report  46  (Feb  1961). 


1 


4 

More  recently,  Greenway  and  Wood  (1973)  performed  a series  of  flow- 
visualization  experiments  to  establish  the  flow  structures  at  the  trailing 
edges  of  stationary  and  oscillating  flat  struts.  In  those  measurements 
the  edges  were  blunt  and  unsymmetrlcally  beveled  with  various  Included 
angles.  The  effect  of  beveling  was  found  to  destabllze  the  vortex 
structure  In  the  wake  but  not  to  Influence  greatly  the  strengths  of  shed 
vortices.  Greenway  and  Wood  calculated  strengths  of  shed  vortices  from 
particle  trajectories  In  photographs.  For  Included  angles  of  less  than 
30  deg,  regularity  In  the  wake  was  found  to  be  nearly  absent  for  rigid 
struts.  These  visualization  results  corroborated  with  the  Heskestad  and 
Olberts  hydroelastic  data  and  led  to  speculation  by  Greenway  and  Wood 
that  a reduction  In  vibration  amplitude  arose  because  of  randomization  of 
the  periodic  forces  on  the  strut  rather  than  from  a reduction  in  the 
mean-square  amplitude  of  the  forces. 

The  Influence  of  splitter  plates  In  altering  the  formation  of 
periodic  vortex  streets  In  the  wakes  of  flat  struts  was  examined  by 
Bearman  (1965).^  Splitter  plates  Interfered  with  Interaction  of  separated 
free-shear  layers  at  bluff  trailing  edges.  Measurements  of  velocity 
fluctuations  in  the  wakes  for  increasing  splitter-plate  lengths  demon- 
strated that  the  disturbances  reached  maxima  at  Increasing  distances 
downstream  from  the  edge.  Locations  of  maximal  disturbances  were  found 
by  smoke  visualization  to  coincide  with  points  between  the  positions  of 
the  first  and  second  vortices.  For  long  enough  splitter  plates, 

Bearman^  found  that  discrete  vortex  formation  could  be  eliminated. 

The  frequency  of  vortex  formation  in  wakes  formed  by  blunt  trailing 
edges  was  first  examined  by  Gongwer  (1952),^  who  proposed  a particularly 
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on  Vortex  Shedding  and  Vibration,"  Journal  of  Fluid  Mechanics,  Vol.  61, 
pp.  323-335  (1973). 

^Bearman,  P.W.,  "Investigation  of  the  Flow  Behind  a Two-Dimensional 
Model  with  a Blunt  Trailing  Edge  and  Fitted  with  Splitter  Plates,"  Journal 
of  Fluid  Mechanics,  Vol.  21,  pp.  2A1-255  (1965). 

^Gongwer,  G.A.,  "A  Study  of  Vanes  Singing  In  Water,"  Journal  of  Applied 
Mechanics,  Vol.  19  (1952). 
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simple  form  of  dimensionless  frequency  or  Strouhal  number.  He  found  by 
experiment  that  the  shedding  frequency  scaled  well  on  the  free-stream 
velocity  and  a length  that  was  the  sum  of  the  edge  thickness  and  the 
momentum  thicknesses  of  the  boundary  layers  forming  the  wake.  Later, 
Bearman  found  for  a blunt  edge  with  a given  length  of  splitter-plate  that 
the  Strouhal  number  formed  with  the  trailing-edge  thickness  and  free- 
stream  velocity  was  nearly  constant  over  a wide  speed  range.  He  found, 
however,  that  the  Strouhal  number  varied  considerably  as  the  splitter- 
plate  length  was  increased.  He  also  found  that  the  static  base-pressure 
coefficient  varied  with  splitter  length  but  not  with  Reynolds  number.  A 
relationship  between  shedding  frequency  and  static  base  pressure  was 
shown  by  Bearman  (1967).^  In  that  paper  Bearman  noted  that  earlier 

O 

definition  of  Strouhal  number  by  Roshko  (1954)  did  not  fit  his  1965 

8 

results  without  modification.  The  definition  by  Roshko  was  based  on  the 
estimated  distance  between  parallel,  separated,  free-shear  layers  in  the 
wake.  This  distance  was  analytically  determined  by  the  notched  hodograph 
method  from  the  potential  flow  around  bluff,  two-dimensional  cylinders. 

The  modification  by  Bearman  replaced  this  distance  with  the  estimated 
cross  wake  separation  between  vortex  centers  in  the  fully  developed  far 
wake.  This  vortex  separation  was  determined  by  use  of  the  Kronauer 

9 

(1964)  hypothesis  that  a minimum  vortex  drag  coefficient  is  exerted  on 
the  wake-producing  body  and  this  coefficient  uniquely  determines  the 
spatial  arrangement  of  vortices.  From  measured  streamwise  vortex  sep- 
arations, Bearman^  calculated  the  cross  wake  separations.  His  defined 
Strouhal  number  appeared  to  be  a universal  function  of  the  mean  base- 
pressure  coefficient  for  data  collected  by  a number  of  investigators  on 


^Bearman,  P.W.,  "On  Vortex  Street  Wakes,"  Journal  of  Fluid  Mechanics, 
Vol.  28,  pp.  625-641  (1967). 
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Kronauer,  R.E.,  "Predicting  Eddy  Frequency  in  Separated  Wakes,"  Paper 
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Symposium  on  Concentrated  Vortex  Motions  in  Fluids,  University  of  Michigan, 
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various  blunt  bodies  with  and  without  splitter  plates.  In  the  definitions 

8 7 

made  by  both  Roshko  and  Bearman,  the  characteristic  velocity  Is  that 
existing  just  outside  the  boundary  layers  on  the  body  at  the  point  of 
separation. 

These  Investigations  have  shown  how  details  of  the  wake  structure 
determine  vortex  shedding  frequencies  and  how  tralllng-edge  shape  In- 
fluences the  Interactions  of  separated  shear  layers.  V/hen  the  shear 
layers  are  decoupled,  vortex  formation  appears  to  become  disordered;  this 
disorder  modifies  the  frequency  of  formation  of  vortices  as  well  as  the 
base  pressure  or  the  drag  coefficient  of  the  resultant  vortex  street. 

What  these  Investigations  do  not  disclose,  however.  Is  the  level  and 
time  dependence  of  dynamic  pressures  which  are  exerted  on  a shedding 
body  for  the  various  tralllng-edge  shapes.  Many  Investigations  of  the 
oscillatory  lift  coefficients  on  shedding  cylinders  have  been  made,  but 
very  little  Information  exists  for  streamlined  aerofoils  with  vortex 
shedding.  Recently,  however,  Hanson  (1970)^^  measured  the  mean  and 
fluctuating  base  pressures  as  well  as  the  near-wake  characteristics  behind 
lifting  airfoils.  His  results  demonstrated  that  a relationship  exists 
among  the  mean  base  pressure  coefficient,  the  distance  of  vortex  formation 
downstream  of  the  edge  SL^,  and  the  momentum  thickness  of  the  wake  6.  The 
mean  base  pressure  coefficient  was  found  to  Increase  linearly  with  the 
ratio  0/il^.  Apparently,  pressure  recovery  of  the  flow  near  the  edge  Is 
associated  with  an  Increase  In  the  formation  length.  Hanson  also 
found  that  the  vortex  shedding  frequency  could  be  scaled  on  momentum 
th:':ckness  of  the  wake  more  adequately  than  a geometric  length,  character- 
istic of  the  trailing  edge. 

This  report  Is  a description  of  measurements  of  fluctuating  pressures 
on  trailing  edges  of  flat  struts  and  of  fluctuating  velocities  In  wakes 
that  generate  those  pressures.  The  measurements  were  made  on  a family 
of  edge  shapes,  similar  to  those  of  Heskestad  and  Olberts^  and  of 


Hanson,  C.E.,  "An  Investigation  of  the  Near-Wake  Properties  Associated 
with  Periodic  Vortex  Shedding  from  Airfoils,"  Massachusetts  Institute  of 
Technology  Acoustics  and  Vibration  Laboratory  Report  76234-5  (1970). 


4 


Greenway  and  Wood.  The  edges  (Figure  1)  were  selected  for  the  range  of 
random  to  periodic  wakes  that  they  produced.  Cross  spectral  densities  of 
the  pressures  on  the  edges  as  well  as  cross  correlations  and  cross 
spectra  of  those  pressures  with  velocity  disturbances  In  the  shear  layers 
near  the  edges  will  demonstrate  the  relationship  of  the  pressures  with  the 
near  wake-formation  regions  of  each  edge.  Chordwlse  distributions  of 
fluctuating  pressure  will  be  shown,  and  Strouhal  number  definitions  will 
be  reviewed.  For  the  trailing  edges  considered,  interrelationships  will 
be  drawn  among  those  which  generate  periodic  surface  pressures  and  those 
which  do  not. 


2.  FACILITY  AND  INSTRUMENT  DESCRIPTIONS 

All  measurements  were  performed  in  the  Anechoic  Flow  Facility  (AFF) 

of  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center  (the 

Center) . This  facility  combines  a low-turbulence  wind  tunnel  with  an 

anechoic  chamber.  The  broadband  turbulence  Intensity  level  in  the  tunnel 

at  a centerline  speed  of  150  ft/sec  has  been  shown  by  DeMetz  and  Casarella 

(1973)^^  to  be  about  0.08  percent.  The  test  section  Is  21  feet  and 

1 Inch  in  length  and  has  a cross-section  of  8 by  8 feet  with  corner  fillets 

to  maintain  one-dimensional  flow  throughout.  Pressure  gradients  which 

are  caused  by  turbulent  boundary-layer  growth  are  reduced  by  tapering  the 

test  section  walls.  Flow  speeds  may  be  controlled  within  +0.5  percent 

with  a maximum  velocity  of  200  feet  per  second.  Measurements  have  been 

made  at  the  end  of  the  test  section  where  It  Is  jointed  by  an  anechoic 

chamber  designed  for  a 150-Hz  lower  acoustic  cutoff  frequency.  A 

description  of  the  acoustic  and  aerodynamic  characteristics  of  the 

12 

facility  has  been  given  by  Bowers  (1973). 
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2 . 1 STRUT  AND  TRAILING  EDGES 

The  working  strut  (Figure  1)  had  a circular  leading  edge  which,  due 
to  laminar  separation,  tripped  the  turbulent  boundary  layer  on  the 
downstream  flat  section.  Excluding  the  removable  tralllng-edge  sections, 
the  strut  was  3 feet  In  chord  with  a unlfon.!  thickness  of  2 Inches  and  a 
span  of  A feet.  Four  trailing  edges  were  used;  on3  was  squared  off,  and 
three  were  unsymmetrlcally  beveled.  Of  the  beveled  edges,  two  had  an 
Included  tip  angle  of  25  deg  and  a length  of  6 1/4  Inches;  the  other  had 
an  Included  tip  angle  of  45  deg  and  a length  of  4 Inches.  The  45-deg 
edge  and  one  of  the  25-deg  edges  were  faired  with  5-  and  10- Inch  radii, 
respectively,  to  give  a continuously  increasing  pressure  gradient.  The 
other  25-deg  edge  had  a knuckle,  2-lnches  from  the  strut  attachment 
point;  see  Figure  1.  Both  the  strut  and  the  trailing  edges  were  con- 
structed from  pine  and  plywood  and  were  painted  to  maintain  a smooth 
surface. 

The  strut  was  rigidly  fastened  In  a vertical  position  to  the  floor 
of  the  tunnel  at  the  end  of  the  upstream  test  section,  adjacent  to  the 
anecholc  chamber.  An  end  plate  In  the  form  of  a plywood  panel  8 feet 
long  and  5 feet  wide  was  attached  to  the  top  of  the  strut  and  to  the 
walls  of  the  tunnel.  The  end  plate  and  tunnel  floor  together  with  the 
vertical  strut  resembled  a sidewise  H.  The  arrangement  effectively 
divided  the  tunnel  In  half.  The  measured  values  of  fluctuating  wall 
pressures  on  the  blunt  trailing  edge  of  the  strut  were  uniform  over  an 
8-lnch  spanwlse  section,  showing  that  the  flow  was  essentially  two 
dimensional  across  the  span  of  the  strut. 

2 . 2 INSTRUMENTATION 
Static  Pressure  Measurements 

The  static  pressure  variation  along  the  chord  was  measured  using 
pressure  taps  at  discrete  points  along  both  sides  of  the  strut.  Including 
the  stagnation  point.  Depending  on  the  trailing  edge,  there  were  up 
to  22  tap  locations.  Each  pressure  tap  consisted  of  1-lnch  long 
and  1/8-inch-OD,  1/16-inch-ID,  copper  tubing,  glued  into  position  so  that 
Its  tip  was  sealed  flush  with  the  surface  of  the  strut.  Removable  panels 
on  one  side  of  the  strut,  permitted  attachment  of  flexible  plastic  tubes 
which  ran  down  through  the  strut  and  under  the  tunnel  floor  to  a bank  of 


manometers. 
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Surface  Pressure  Fluctuation  Measurements 


The  flow- induced,  surface-pressure  fluctuations  on  the  strut  were 
measured  using  a 1/8-inch  condenser  microphone  (Bruel  and  Kjaer  (B  & K) 
Model  4138)  whose  sensing  area  was  reduced  by  placing  a cap  perforated 
with  a 1/32-inch-diameter  hole  in  its  center  over  the  diaphragm  of  the 
microphone.  The  1/8-inch  condenser  microphone  was  connected  to  the 
cathode  follower  (B  & K Model  2619)  through  a 1/4-lnch,  right-angle 
adaptor  (B  & K Model  UA  0122)  which  was  modified  to  accept  the  1/8-lnch 
microphone.  The  pinhole  perforating  the  center  of  this  cap  led  to  a 
small  cavity  above  the  diaphragm.  The  frequency  response  of  the  pinhole 
microphone  was  measured  by  comparing  the  response  of  the  microphone  with 
and  without  pinhole  caps  to  identical  free-field  sound  fields.  The 
relative  phase  shift  between  the  microphones  was  also  measured  in  this 
manner.  The  Helmholtz  resonance  frequency  caused  by  the  cavity  was 
measured  at  17.5  kilohertz.  Response  of  the  pinhole  microphone  was  from 
flat  to  1 dB  from  65  Hz  to  10  kHz,  and  the  relative  phase  shift  was  from 
zero  to  5 kHz,  5 deg  at  8 kHz,  and  14  deg  at  10  kilohertz. 

The  1/8-lnch  condenser  microphones  with  pinhole  caps  were  flush 
mounted  to  various  positions  along  the  strut  and  trailing  edges  at  the 
positions  indicated  in  Figure  1.  The  autospectra  as  well  as  the  longi- 
tudinal (streamline)  and  lateral  (spanwise)  cross  spectral  densities  of 
pressures  were  measured  at  these  locations.  The  dimensions  of  the 
microphone  positions  will  be  described  as  the  measurements  are  presented. 

Besides  the  1/8-inch  microphone,  a 1/2-inch  microphone  (B  & K 
Model  4133)  with  a probe-tube  extension  (B  & K Model  UA  0040)  was  used  to 
measure  the  pressure  spectra  at  narrow  thicknesses  of  the  trailing  edge 
as  well  as  the  tip.  The  probe  tube  acts  as  a low-pass  filter  between 
the  microphone  and  the  fluctuating  pressure  source.  Three  3/64-inch- 
dlameter  probe  tubes  were  calibrated  with  lengths  of  1 1/2,  3,  and 
4 7/16  inches;  all  obtained  a fldC  amplitude  response  to  2 kHz  and 
approximately  a linear-phase  response.  Calibration  of  the  probe  tube 
with  proper  damping  was  accomplished  by  using  a B & K coupler  that 
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accepted  a probe  tube,  a 1/2-inch  microphone,  and  a small  sound  source. 

The  volume  of  the  coupler  is  much  smaller  than  the  wavelength  of  sound, 
thereby  producing  a uniform  sound  pressure  level  in  the  chamber  over  the 
frequency  range  of  Interest.  The  source  was  driven  by  a frequency 
analyzer  (General  Radio  Model  1900)  connected  in  the  tracking  mode.  The 
signals  from  both  the  1/2-inch  microphone  and  the  probe  tube  were  fed 
back  into  the  analyser  for  frequency  analysis  and  were  displayed  on  a 
graphic  level  recorder  (General  Radio) . The  response  of  the  tube  was 
adjusted  by  placing  damping  material  (steel  wool)  Into  the  probe  tube 
and  comparing  the  amplitude  and  phase  response  of  the  probe  tube  with  that 
of  the  1/2-inch  microphone.  Responses  are  shown  in  Figures  2 and  3. 

Velocity  Measurements 

Turbulent  velocity  fluctuations  were  measured  using  a heat  flux 
system  (Thermo-Systems  Model  1010),  constant  temperature  anemometer,  and 
a llnearlzer  (Thermo-Systems  Model  1005B).  The  anemometer  was  calibrated 
for  wind  speeds  using  a venturi  meter.  A 0. 0005-Inch,  platlnum-irldlum 
alloy  wire  with  a 0.3-msec  time  constant  was  used  for  the  sensing  element 
which  was  attached  to  an  18-lnch  probe  support.  The  probe  support  was 
positioned  by  an  electromechanical  linear  actuator  to  within  0.01  inch, 
+0.005  inch,  by  reading  a dc  voltage,  fed  back  from  a linear  potentiometer. 
Probe  and  probe  supports  were  aerodynamlcally  shaped  to  minimize  vibration 
from  vortex  shedding  from  the  supports.  Turbulence  intensities  were 
measured  using  a B & K true  rms  meter. 

Auto  and  Cross  Spectrum  and  Correlation 
Analyses 

All  of  the  wall  pressure  auto  and  cross  spectral  densities  as  well 
as  the  pressure-velocity,  cross  spectral  densities  were  obtained  with  a 
real-time  analyzer  system  (Time  Data  Type  1923/A) . The  spectral 
densities  were  determined  digitally,  using  a fast  Fourier-transform 
algorithm  with  a Hanning  window.  Sample  sizes  of  at  least  256  were  used 


with  effective  analysis  bands,  ranging  from  12.5  to  62.5  hertz.  The 
bandwldths  of  analyses  were  constant  over  each  complete  frequency  range, 
which  was  either  from  0 to  2 or  from  0 to  20  kilohertz.  Some  velocity 
spectral  densities  were  obtained  In  1/10-octave  bands  using  a sound  and 
vibration  analyzer  (General  Radio  Model  1564). 

Broadband  correlations  were  made  with  a digital  correlator  (Hewle,.t 
Packard  Model  3721A).  This  Instrument  displayed  both  auto  and  cross 
correlation  functions  at  100  computed  points  with  time  delays  from 
100  psec  to  100  sec  possible.  The  digital  equivalent  of  resistor- 
capacitor  circuit  (RC)  averaging  was  used. 

2.3  STATISTICAL  PRECISION 

Figure  4 shows  measured  and  theoretical  expressions  for  the  standard 
deviation  of  cross  spectral  density  estimates  of  random  signals.  The 
measurement  was  conducted  by  correlating  two  broadband  signals,  which  were 
known  to  be  uncorrelated,  using  the  real-time  analyzer.  The  computed 
values  of  the  normalized  cross  spectrum,  theoretically  zero,  were  estimate 
errors  which  decreased  with  Increasing  sample  size  N.  For  each  sample 
size  a value  of  the  computed  cross-spectrum  was  selected.  This  value 
exceeded  80  percent  of  the  estimates  of  the  spectrum,  and  it  was  In- 
terpreted as  an  80-percent  confidence  limit  on  the  error  of  a cross 
spectral  density  estimate.  A theoretical  estimate,  whose  derivation  Is 
presented  In  Appendix  A,  Is 
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<j)^^(f)  Is  the  true  autospectral  density  of  the  Input 

Is  the  true  cross  spectral  density 

Is  the  standard  deviation  of  the  cross  spectrum 
estimate . 


3.  CHARACTERISTICS  OF  FLOW  INCIDENT  ON  EDGES 
The  measurements  described  in  this  section  served  to  establish  that 
the  turbulent  flow  field,  Incident  on  the  edges,  was  typical  of  that 
which  exists  in  fully  developed  boundary  layers  on  flat  surfaces. 
Measurements  of  the  correlations  between  fluctuating  wall  pressures  at 
different  chordwlse  points  as  well  as  of  turbulent  velocity  were  made. 
These  compared  favorably  to  the  measured  statistics  of  other  Investigators; 
they,  therefore,  confirmed  proper  experimental  technique,  and  standardized 
the  Incident  flow. 

Figure  5 shows  the  chordwlse  static  pressure  distribution  and  the 
cross-section  shape  of  the  strut,  which  was  outfitted  with  the  25-deg — 
rounded  trailing  edge.  The  static  pressures  are  expressed  as  a 
coefficient 
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which,  in  the  absence  of  a local  pressure  gradient  across  the  boundary 
layer,  can  be  interpreted  in  terms  of  the  local  free-stream  velocity  on 
the  strut  U^.  In  Equation  (1)  is  the  ambient  pressure  far  from  the 
strut,  P is  the  local  static  pressure  on  the  strut,  p is  the  fluid 
density,  and  is  the  mean  velocity  upstream  of  and  Incident  on  the 
strut.  Locations  1 through  6 In  Figure  5 designate  single  microphone 
locations  for  which  the  static  pressure  distribution  was  trivially 
altered  by  modifying  the  edges.  Measurements  at  all  these  locations  were 
made  only  once  In  the  course  of  the  experiment,  and  representative  checks 
were  made  for  each  edge.  The  lettered  locations  near  the  trailing  edge 
are  those  for  which  the  static  pressure  was  influenced  to  varying  degrees 
by  the  different  edges;  measurements  were  performed  at  those  locations  for 
all  edges.  The  edges  were  attached  to  the  main  strut  36  Inches  from  the 
leading  edge. 


The  results  in  Figure  5 show  a sharp  static  pressure  minimum  forward 

of  microphone  position  1,  followed  by  a steep  adverse  pressure  gradient. 

Flow  visualization  with  oil  streaks  at  100  ft/sec  disclosed  a separation 

region  extending  from  approximately  1 Inch  to  3 Inches  downstream  from 

the  leading  edge.  It  Is  conjectured  that  laminar  separation  was  caused 

by  the  circular  leading  edge.  Short  separation  bubbles  on  leading  edges 
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of  airfoils  as  discussed  by  Chang  (1970)  play  the  role  of  causing 
transition  for  the  downstream  boundary  layer.  In  this  type  of  flow  the 
separated  laminar  shear  layer  becomes  turbulent  before  reattachment. 
Likely,  the  free-stream  turbulence  level  of  the  wind  tunnel  (0.08  per- 
cent) was  low  enough  to  permit  laminar  separation  to  occur.  The  flow 
around  the  cylindrical  leading  edge  In  the  present  case  Is  probably  some- 
what similar  to  that  around  the  forward  quadrant  of  a circular  cylinder. 

13 

In  that  case  (Chang)  laminar  separation  occurs,  82  deg  removed  from  the 
forward  stagnation  point  for  Reynolds  numbers,  based-diameter,  less  than 
2 X 10^.  Turbulent  separation  occurs  further  around  for  Reynolds  numbers 
greater  than  7 x 10^.  In  the  present  Instance,  for  * 100  ft/sec  the 
Reynolds  number  based  on  the  diameter  of  the  leading  edge,  h > 2 Inches, 
Is  1.02  X 10^.  The  conjecture  of  laminar  separation,  then,  seems 
reasonable. 

Further  downstream,  the  static  pressure  varies  slightly,  becoming 
somewhat  favorable  forward  of  the  trailing  edge.  Magnitudes  of  the 
favorable  and  adverse  pressure  gradients  at  the  trailing  edge  are 
dependent  on  the  shape  of  the  trailing  edge;  each  edge  will  be  discussed 
In  detail  later  In  the  report.  For  all  except  the  blunt  trailing  edge, 
the  favorable  pressure  gradient  from  x > 28  Inches  to  x « 34  Inches 
shown  In  Figure  5 Is  typical.  Therefore,  the  statistics  of  the  boundary 
layer  for  all  positions  1 through  B are  expected  to  be  unaltered  by 
changing  the  beveled  trailing  edges  but,  perhaps,  slightly  modified  by 
the  blunt  edge.  This  difference  In  local  static  pressure  Is  thought  to 
be  caused  by  the  slight  change  In  camber  brought  about  by  the  use  of 
unsymmetrlcal  versus  symmetrical  edges. 
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Chang,  P.K. , "Separation  of  Flow,"  Permagon  Press,  Oxford,  England 


The  mean  and  turbulent  velocity  profiles  for  positions  ^and  B are 

shown  in  Figure  6.  The  mean  velocity  U(y)  and  root-mean-square  turbulent 

-jl/2 

velocity  u are  expressed  as  fractions  of  the  local  free-stream 

velocity  U^.  Distance  normal  to  the  plane  of  the  strut  y is  normalized 

on  the  boundary- layer  thickness  6,  which  is  defined  by  U(6)  = 0.99  U^. 

The  mean  velocity  profiles  appear  to  be  similar,  indicating  that  an 

equilibrium  boundary  layer  has  been  formed  on  the  strut.  The  Prandtl 

14 

(Schlichting  (1960)  ) 1/7-power  law  for  the  mean  velocity  distribution 

is  also  shown,  and  it  compares  favorably  with  the  results.  This  compari- 
son further  shows  that  the  mean  velocity  in  the  boundary  layer  is  similar 
to  that  on  a flat  plate.  Also,  the  Klebanoff  (1955)^^  longitudinal 
turbulence  intensities  are  shown;  levels  in  the  outer  region  of  the 
boundary  layer  on  the  strut  are  higher  than  those  of  Klebanoff.  That 
they  decrease  from  position  4 to  position  B suggests  that  these  levels 
are  caused  by  large-scale  disturbances  associated  with  the  leading-edge 
separation.  Close  to  the  wall  y/6  < 0.15,  the  measured  intensities  are 
lower  than  those  measured  by  Klebanoff. The  discrepancy  could  be  due 
to  a probe-wall  interference  or  possibly  (but  less  likely)  to  the  effect 
of  upstream  history  determined  by  the  leading-edge  separation. 

Spectral  densities  of  the  turbulent  velocities  at  selected  locations 
in  the  boundary  layer  at  position  B on  the  knuckle  edge  are  shown  in 
Figure  7.  The  spectra  are  normalized  on  the  local  mean  velocity  U(y), 
the  boundary-layer-displacement  thickness 


6 


and  the  mean  square  velocity  u . They  are  two  sided  so  that 
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Schlichting,  H. , "Boundary  Layer  Theory,"  McGraw-Hill  Fourth  Edition, 
New  York  (1960) . 

^^Klebanoff,  P.S.,  "Characteristics  of  Turbulence  in  a Boundary  Layer 
with  Zero  Pressure  Gradient,"  National  Advisory  Committee  for  Aeronautics 
Report  1247  (1955). 
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(3) 


■ J 


r(i)6*  1 

Lu(y)J 

where  w is  the  radian  frequency.  Near  the  wall  for  y/6  = 0.028,  the 

normalized  spectrum  is  deficient  in  low-frequency  intensity,  compared  to 

spectra  further  from  the  wall.  Generally  the  spectra  are  comparable  to 

15  16 

those  measured  elsewhere,  e.g.,  Klebanoff,  and  Blake  (1970). 

Tables  1 and  2 give  mean  boundary-layer  properties,  measured  else- 
where on  the  strut  with  25-deg  knuckle  and  round  trailing  edges.  The 
momentum  thicknesses 


and  the  'jhape  factors 

H = 6*/e  (5) 

are  shown  with  6 and  6*  for  each  location.  At  position  B for  100  feet/sec 
the  boundary  layer  on  the  25-deg  knuckle  edge  is  slightly  thicker  than 
that  at  the  same  location  for  the  25-deg  round  edge. 

The  mean  velocity  profiles  at  positions  A and  B forward  of  the  25- 
deg  round  trailing  edge  are  plotted  logarithmically  for  100  ft/sec  in 
Figure  8.  Velocities  are  normalized  on  the  friction  velocity 


(6) 

where  T is  the  local  wall  shear.  The  value  of  U was  determined  by  curve 

W T ^ 

I - - 

fitting  to  the  law  of  the  wall  as  described  by  Perry  and  Joubert  (1963). 


u 


I 


/ 


Blake,  W.K.,  "Turbulent  Boundary  Layer  Wall  Pressure  Fluctuations  on 
Smooth  and  Rough  Walls,"  Massachusetts  Institute  of  Technology  Acoustics 
and  Vibration  Laboratory  Report  70208-1  (1969). 

^^Perry,  A.E.  and  P.N.  Joubert,  Journal  of  Fluid  Mechanics,  Vol.  17, 
p.  193  (1963). 
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so  Feet  per  Second 


100  Feet  per  Second 


Pos 

6 

B 

e 

H 

6 

6* 

0 

H 

A 

0.073 

1.23 

B 

0.103 

0.074 

1.4 

C 

0.0705 

0.0564 

1.245 

D 

0.90 

0.084 

0.0795 

1.055 

0.05 

0.0475 

1.06 

E 

0.07 

0.0525 

1.34 

E-1 

1.33 

wmm 

3.4 

0.88 

0.316 

0.114 

2.765 

F 

1.6 

m 

4.47 

1.29 

0.596 

0.217 

2.75 

All  dimensions  are  in  inches,  I 

1 - 6*/e. 

TABLE  2 - MEAN  BOUNDARY-LAYER  PROPERTIES  FOR  25-DEGREE  ROUND  TRAILING  EDGE 


60  Feet  per  Second 

100  Feet  per  Second 

Pos 

6 

6* 

0 

H 

5 

6* 

0 

H 

4 

0.575 

0.0735 

0.061 

1.21 

0.435 

0.0603 

0.048 

1.26 

A 

B 

0.8 

0.093 

0.08 

1.17 

0.63 

0.073 

0.059 

1.23 

B-1 

0.615 

0.0668 

0.0565 

1.16 

C 

0.5 

0.045 

0.04 

1.14 

DE 

0.8 

0.0915 

0.059 

1.56 

0.62 

0.093 

0. 0495 

1.88 

E-1 

0.95 

0.157 

0.082 

1.93 

0.65 

0.118 

0.0676 

1.74 

F 

1.35 

0.285 

0.133 

2.12 

0.95 

0.231 

0.095 

2.43 

G 

1.96 

0.52 

0.215 

2.42 

1.36 

0.451 

0.129 

3.47 

All  dimensions  are  in  inches,  H - 6*/6. 


100  ft/sec 


The  local  wall-shear  coefficient  at  x = 32  inches  and  U„  * 

00 

C[  - \n/2  p (7) 

-3 

was  found  to  be  C,  « 3.2  x 10  , somewhat  lower  than  that  given  by 

14  ^ 

Schlichting  for  a naturally  developed  boundary  layer  on  a smooth  wall 
at  this  Reynolds  number  U^x/v. 

i The  mean-square-wall-pressure  fluctuations,  measured  in  62.5-Hz  fre- 

i "2 

quency  bands  at  positions  1 through  6,  p (f ) , are  shown  in  Figure  9 as 
fractions  of  free-stream  dynamic  head  q^  and  a function  of  frequency  f for 
* 50  feet  per  second.  At  position  1 the  narrowband  pressures  at  low 
frequencies  decrease  from  an  undetermined  maximum  at  f < 120  Hz  and  then 
show  a broad  peak  at  2700  hertz.  The  levels  at  position  2 are  broadband; 
at  locations  which  are  further  removed  from  the  leading  edge,  the  levels 
continue  to  decrease  slightly.  At  100  ft/sec  the  narrowband  wall 
pressures  (Figure  10)  at  position  1 are  more  broadband  than  at  50  ft/sec; 

however,  they  are  substantially  higher  than  at  other  locations  on  the 

_2 

strut.  There  is  a broad  "hump"  centered  on  4500  Hz,  followed  by  f be- 
havior. The  oscilloscope  trace  of  the  unfiltered  signal  showed  that  be- 
tween * 75  and  100  ft/sec,  as  speed  was  increased  slowly  and  continuously, 
I the  pressure  fluctuations  "jumped"  markedly.  This  suggests  that  the 

S disturbances  at  position  1 which  are  associated  with  the  leading-edge 

separation  changed  in  character  as  the  speed  was  increased  from  50  to 
I 100  feet  per  second.  Perhaps  the  separated  shear- layer  reattachment 

occurred  immediately  above  position  1 at  100  ft/sec  and  a little  downstream 
from  position  1 at  50  feet  per  second.  The  presence  of  the  peaks  at 
2700  and  4500  Hz  indicates  the  occurrence  of  pulsation  associated  with 
the  laminar  separation  reattachment  whose  frequency  increases  nearly 
linearly  with  U^.  The  pressure  spectra  at  locations  downstream  from 
position  1 are  more  uniformly  behaved  through  position  B. 

To  determine  how  the  wall  pressures  behave  in  relation  to  earlier 
work,  the  spectra  at  position  B for  2 speeds  have  been  normalized  on  the 
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I 


local  wall  shear  coefficient. 
Figure  11.  Normalization  has 


the  displacement  thickness,  and  Inflow  In 
been  made  so  that 


4>({i))U 


6* 

w 


d 


(8) 
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The  results  compare  favorably  to  those  of  Blake  (1970)  ’ for  a boundary 

layer  on  a smooth,  flat  wall.  For  > 4 the  results  diverge  because 

of  the  spatial  averaging  of  small-scale  disturbances  over  the  microphone 

19 

face  as  discussed  by  Corcos  (1963).  Those  effects  become  most  severe 
when  U)d/U^  > 2,  where  d Is  the  diameter  of  the  microphone,  and  - 0.6 
Is  the  convection  velocity  of  disturbances  across  the  microphone 
considered  In  Section  6.  For  the  current  case,  d > 0.3  5*,  U - 0.6  U 

' C 00 

(Section  4),  and  a)6*/U  ==  4 so  that  o)d/U  ==  2.  The  spectra  also  diverge 

00  Q 

slightly  at  low  dimensionless  frequencies.  This  divergence  was  not 

16  18 

observed  In  earlier  work  (Blake  * ),  and  In  this  case  It  may  be  due  to 

the  Influence  of  the  leading-edge  separation.  It  could  also  be  associated 
with  the  fact  that  the  current  measurements  were  obtained  at  a lower 
Reynolds  number  6/v  than  In  the  earlier  work.  Nonetheless 

®0 

these  dimensionless  spectra  are  considered  to  be  In  good  agreement  with 
the  earlier  ones. 

The  spatial  cross  spectral  densities  of  wall  pressures  at  position 
B were  measured  with  a pair  of  microphones  separated  In  the  streamline 
direction  by  a variable  distance  r^.  The  magnitudes  of  the  cross  spectrum 
|<t>(r^.O,o))|  are  normalized  on  the  autospectrum  <I>(ro)  and  are  plotted 
against  the  phase  a(r  ,0))  In  Figure  12.  For  separations  of  r ■ 36*  and 

X X 

r^  = 65*,  the  magnitude  Is  a function  of  phase;  this  behavior  Is  typical 

of  pressure  fluctuations  on  smooth  walls.  The  dotted  curve  In  Figure  12 

18 

Is  from  the  data  of  Blake  which  were  obtained  on  the  smooth  wall 


18 

Blake,  W.K. , "Turbulent  Boundary-Layer  Wall  Pressure  Fluctuations  on 
Smooth  and  Rough  Walls,"  Journal  of  Fluid  Mechanics,  Vol.  44,  p.  1 (1970). 

19 

Corcos,  G.M.,  "On  the  Resolution  of  Pressure  In  Turbulence,"  Journal 
of  the  Acoustical  Society  of  America,  Vol.  35,  p.  192  (1963).. 
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beneath  a boundary  layer  of  substantially  larger  Reynolds  number.  The 

19 

phase  of  the  cross  spectrum  has  been  Interpreted  by  Corcos  in  terms  of 
a convection  velocity  (note  r^  * r^) 

a(r^,o))  - (9a: 

or  in  terms  of  a characteristic  wavelength 

a(r^,o))  - 2Tir^/X^  (9b] 

Using  this  interpretation,  the  data  in  Figure  12  suggest  that  in  the 
present  case  the  coherence  of  disturbances  diminishes  to  0.1  within  two 
wavelengths,  compared  to  three  wavelengths  beneath  the  boundary  layer  at 
high  Reynolds  number.  The  convection  velocities  which  were  determined 
from  the  phases  will  be  discussed  in  Section  6.  Cross  spectral  densities 
for  a lateral  separation  of  r^  » 36*  are  shown  in  Figure  13.  For 
tiir^/U^  >1.5  they  are  similar  to  those  measured  on  other  smooth  walls; 
for  wr  /U  < 1.5  they  are  less  than  those  measured  by  Blake;^®  however, 
the  discrepanv'.y  may  be  related  to  an  additional  separation  dependence  as 
noted  by  Bull  (1967).^^ 

4.  PRESSURES  AND  VELOCITIES  GENERATED  BY  PERIODIC  VORTEX  STREETS 
Disturbances  in  the  wakes  of  both  the  blunt  trailing  edge  and  the 
rounded  45-deg  trailing  edge  were  characterized  by  the  development  of 
periodic  vortex  streets.  In  this  section  we  will  consider  first  the 
particular  details  of  the  flow  off  each  edge  separately  and  second  the 
general  characteristics  of  the  periodic  disturbances. 

4.1  WAKES  OF  BLUNT  TRAILING  EDGES 

The  wake  structures  behind  blunt  trailing  edges  have  been  determined 


in  detail  by  Bearman^’^  and  Greenway  and  Wood.^ 


Current  measurements  on 


the  blunt  edge  corroborate  and  extend  those  results,  and  they  serve  as  a 
basis  for  interpreting  the  other  measurements  to  be  described  later. 


Bull,  M.K.,  "Wall  Pressure  Fluctuations  Associated  with  Subsonic 
Turbulent  Boundary  Layer  Flow,"  Journal  of  Fluid  Mechanics,  Vol.  28,  p.  719 
(1967). 


The  characteristics  of  the  flow  field  of  the  blunt  trailing  edge  are 
shown  In  Figure  14.  The  descriptions  begin  29  Inches  downstream  from  the 
leading  edge  and  extend  41/2  Inches  downstream  from  the  trailing  edge. 

The  static  pressure  distribution,  shown  at  the  top  of  the  figure,  be- 
comes favorable  near  the  edge;  the  base  pressure  coefficient  which  was 
measured  at  numerous  locations  on  the  base  Including  the  centerline  of 
the  strut  was 

- C = 0.51 
Pb 

It  was  uniform  across  the  thickness  of  the  edge  and  the  span.  The 

Reynolds  number,  based  on  chord  U^w/v,  where  w Is  the  chord,  extended 

to  2 X 10^.  For  the  experiments  of  Bearman^  It  extended  to  2.85  x 10^, 

4 4 

and  In  the  case  of  Greenway  and  Wood  It  extended  to  1.4  x 10  for  their 
flow  visualizations  and  to  1.3  x 10^  for  their  wind  tunnel  measurements. 

The  vortex  shedding  frequency  for  the  current  strut,  measured  both 
with  fluctuating  pressures  on  the  strut  and  with  the  velocity  fluctuations 
In  the  wake,  was  found  to  be  linearly  dependent  on  U^.  To  compare 
shedding  frequencies  In  this  and  earlier  work,  we  note  the  Strouhal 
number  based  on  the  strut  thickness  h,  which  Is 

N ^ - f h/U„  - 0.23  (10) 

St  s <» 

4 

This  number*  Is  somewhat  smaller  than  that  of  Greenway  and  Wood  at  a low 

Reynolds  number  » 0.286;  however.  It  Is  In  agreement  with  their  value 

of  * 0.24  obtained  at  a high  Reynolds  number.  This  difference  Is 

probably  caused  by  laminar  separation  at  the  trailing  edges  at  their  low 

Reynolds  number;  also,  the  shed  boundary  layers  In  this  and  the  Greenway 

and  Wood  cases  are  dissimilar.  The  current  value  compares  favorably  with 

the  Bearman^  value,  N = 0.24,  for  a blunt  edge  without  a splitter  plate. 

s c 

All  of  the  previously  mentioned  differences  In  dimensionless  shedding 


In  this  case  the  definition  Is  practically  equivalent  to  that  of 

Gongwer  which  replaces  h by  h plus  the  sum  of  the  momentum  thicknesses  of 
the  boundary  layers  shed  off  the  edge.  Thus,  h+26*2.1  Inch  and  f (h+20)/ 

= 0.25  which  Is  somewhat  greater  than  the  value  0.19  obtained 

by  Gongwer. 
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frequency  are  most  likely  due  to  variations  In  the  wake  thicknesses 

occurring  In  the  various  experiments,  and  they  reflect  the  observation 

(Bearman)^  that  Equation  (10)  Is  not  a universal  scaling  number.  One  of 

the  purposes  of  this  report  is  to  review  various  frequency-scaling 

hypotheses  In  light  of  the  current  measurements. 

The  mean  velocity  profiles  of  the  flow  on  the  strut  and  in  the  near 

wake  of  the  strut  are  shown  at  the  center  of  Figure  14.  At  the  bottom 

of  Figure  14  are  shown  the  vertical  profiles  of  root-mean-square  filtered 

— 1/2 

fluctuating  velocity  u turbulence  velocities  were  filtered 

Ip  1/10-octave  bands  centered  on  the  shedding  frequency  f . Within 

s 

0.3  inch  of  the  strut  the  profiles  display  sharp  maxima  near  the  corners 

of  the  edge  and  broader  peaks  closer  to  the  centerline  of  the  wake.  The 

sharp  peaks  In  intensity  are  associated  with  the  primary  disturbances  In 

the  shear  layers,  while  the  peaks  near  the  centerline  are  probably 

caused  by  a secondary  vortex  system.  For  x > 0.5  inch,  the  profiles  of 

— 1/2 

u (f^)  display  two  maxima  which  broaden  In  spatial  extent  as  the  dis- 
turbances are  convected  away  from  the  strut.  The  locus  of  the  filtered 
Intensity  maxima  Is  sketched  In  the  center  of  Figure  14;  It  Is  situated  a 
distance  which  we  shall  call  y « above  the  centerline  of  the  wake. 

The  root-mean-square  fluctuating  pressure,  filtered  In  a 12.5-Hz 

“2 

band  centered  on  the  shedding  frequency  p (f  ) > was  measured  at 

s 

increasing  distances  from  the  trailing  edge.  The  pressures,  normalized 

on  the  dynamic  head  of  the  Inflow,  are  shown  In  the  upper  part  of 

Figure  14.  The  fluctuating  pressures  on  one  side  of  the  strut  were  In 

common  phase,  showing  no  convection  as  was  expected.  The  frequency 

spectra  of  fluctuating  pressure  on  the  surface  of  the  strut  (Figure  15) 

displayed  a single  dominant  peak  at  f > f^.  (This  figure  will  be  more 

fully  described  later.)  An  autospectrum  of  the  fluctuating  base  pressure 

displayed  a single  peak  at  f ■ 2f  . The  narrowband,  root-mean-square 

s 

base  pressure,  filtered  at  f ■■  2f^,  is  also  shown  In  Figure  14. 


19 


The  growth  of  coherent  disturbances  In  the  near  wake  can  be  deduced 
from  correlations  of  fluctuating  pressure  on  the  edge  with  fluctuating 
velocity  at  specific  locations  In  the  wake.  The  magnitudes  of  the 
fluctuating  pressures  generated  by  the  wake  should  depend  on  both  the 
strength  of  the  vortices  shed  into  the  wake  and  the  proximity  of  the 
vortices  to  the  edge.  The  frequency  of  the  pressures  is  the  same  as  the 
rate  of  formation  of  vortices,  and  this  Is  given  by  the  ratio  of  the  local 
speed  of  convected  vortices  relative  to  the  fixed  edge  to  the  streamwlse 
spacing  of  adjacent  vortices  in  a single  row.  The  cross-wake  spatial 
variation  of  local  fluctuating  velocity  In  a vortex  street  has  been  related 
to  the  geometry  of  the  configuration  of  the  vortices  by  Schaefer  and 

21  ~~2 
Esklnazl  (1959).  They  have  shown  that  the  maximum  of  u (f  )>  l.e., 

s 

y ■ with  the  locus  of  the  peripheries  of  the  vortex  cores 

rather  than  the  locus  of  centers  of  the  vortices.  The  diameter  of  each 
vortex  core  was  shown  to  extend  approximately  from  the  position  of 
maximum  filtered  velocity  to  the  centerline  of  the  wake.  The  strength  of 
the  vortex  Is  proportional  to  the  product  of  the  diameter  of  the  vortex 
core  and  the  tangential  velocity  at  the  extremity  of  the  core.  Therefore, 
measurements  of  the  correlation  between  the  fluctuating  surface  pressure 
and  the  fluctuating  velocity  along  the  y = are  expected  to  reflect 

the  structure  of  vortles  In  the  wake. 

This  measurement  was  performed  by  moving  an  anemometer  probe  along 
the  locus  of  filtered  turbulence  maxima  In  the  upper  side  of  the  wake 
that  Is  shown  In  Figure  14.  The  space-time  correlation  of  pressure  and 
velocity  is  written 


<r,.T)  ■ 


-0.375,  t)  u(r^,i:  + T) 


U 

o 


(11) 
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Schaefer,  J.W.  and  S.  Esklnazl,  "An  Analysis  of  the  Vortex  Street 
Generated  in  a Viscous  Fluid,"  Journal  of  Fluid  Mechanics,  Vol.  6,  p.  241 
(1959). 
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where  p (t)  “ p(f  ■ -0.375  ln.,t)  is  the  unflltered  fluctuating  pressure 

as  close  as  possible  to  the  tip 

r “ -0.375  inch 
^1 

t Is  time 

u(r^,t)  is  the  local  velocity  along  y = y^l'2- 

r^  Is  measured  downstream  of  the  trailing  edge. 

The  correlation  has  been  normalized  on  the  root-mean-square  fluctuating 

pressure  at  the  tip  and  the  local  free-stream  velocity  U . By  using  U 

~2  1/2  ° ° 

rather  than  the  local  root-mean-square  velocity  u (f  ) as  a normalizing 

• 8 

factor,  we  can  determine  the  Influences  of  specific  locations  of  the  wake 
on  the  fluctuating  pressure.  The  measured  correlations  for  = 100  ft/sec 
are  shown  In  Figure  16.  The  upper  curve  Is  the  spatial  correlation,  which 
reaches  a value  R (r  ,o)  = -0.25  at  r =0.5  inch  and  a value  R (r  ,o)  = 

pu  X X pu  X 

+0.25  at  r =2.1  Inch.  Zero  crossings  occur  at  r =1.1  inch  and 
* ^ 1/2  ^ 
r =3.5  inch.  The  maximum  values  of  u (f  )/U  , taken  from  Figure  14, 

^ ° ~2  1/2 

show  an  absolute  maximum  at  r =1.3  Inch;  the  intensities  grow  from  u 

^ 5 

■ 0.16  U to  a maximum  0.29  U . This,  according  to  the  Bearman  criterion, 
o 0 

signifies  the  end  of  the  vortex  formation  region.  The  correlation 

coefficient  evaluated  at  constant  phase  Is  determined  by  R (r  ,T  ),  where 

pu  X m 

T Increases  from  zero,  and  it  essentially  correlates  the  pressure  with  a 
in 

single  vortex  as  it  moves  downstream.  We  shall  call  this  the  moving-eddy 

correlation,  and  It  shows  that  although  the  coherence  is  greatest  at  the 

end  of  the  formation  region  r = its  magnitude  essentially  follows 

— 1/2  X t _ ^^2 

that  of  u (f  )•  The  correlation  coefficient  based  on  u (f  ) > 

s c 

therefore,  is  unity  which  simply  means  that  the  pressure  is  caused  entirely 
by  a fully  coherent  vortex  street.  The  values  of  the  time  delay  for 
the  moving-eddy  correlation  are  shown  at  the  location  of  Figure  17.  The 
varying  slope  in  the  curve  of  X versus  r shows  that  when  the  vortex 
begins  to  form  from  the  free  shear  layer,  the  pressure  is  generated 
Instantaneously  on  the  surface  of  the  strut.  As  the  vortex  continues  to 


"iMirr  II 
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form,  it  begins  to  move  downstream  from  the  strut  while  Its  strength 
Increases.  This  acceleration  continues  until  r *1.3  Inch,  at  which 

X 

point  the  acceleration  diminishes,  and  It  begins  to  move  with  a slowly 

varying  speed  of  convection.  The  convection  velocity,  defined  as 

» Ar  /At  at  this  point  Is  0 = 0.64  U . This  convection  velocity,  based 

X m c o 

on  the  moving-eddy  correlation,  was  also  obtained  by  using  a pair  of  hot 
wire  anemometers  and  forming  velocity  space-time  correlations.  The 
points  obtained  In  the  near  wake  have  been  Included  In  part  d of 
Figure  16,  and  they  fully  corroborate  the  pressure-velocity  correlation. 
In  the  far  wake,  l.e.,  for  ^ Inches  or  ^j^/h  * 4.5,  the  velocity 

correlations  give  0^  * 0.81  U^.  According  to  the  Bearman^  measurements 
of  velocity  correlations,  the  convection  velocity  of  disturbances  In  the 
wake  became  constant  and  equal  to  0.88  U for  r /h  > 4. 

00  X 

To  construct  a geometric  picture  of  the  vortex  structure  we  use  the 
convection  velocities,  determined  from  the  data  of  Figure  16,  and  the 
frequency  of  shedding  f^  to  determine  streamline  vortex  separations  b 
from  the  relationship 

U ■ f b 
c s 

Convection  velocities  are  shown  at  the  bottom  of  Figure  17.  In  the  near 
wake,  the  spatial  correlation  of  Figure  16a  Is  Interpreted  to  show  that 
the  first,  growing  vortex  Is  formed  at  “ 0.5  Inch  duvmstream  from  the 
edge.  At  the  same  time,  the  vortex,  shed  one-half  cycle  earlier  from  the 
opposite  corner  of  the  edge  and  grown  to  full  strength,  has  moved  to 
r * 2.1  Inches.  The  cross  wake  position  of  the  vortex  centers  In  the 

X 21 

near  wake  can  be  estimated  by  using  the  Schaefer  and  Esklnazl  (1959) 
semlemplrlcal  result,  which  roughly  locates  the  vortex  cores  between 
y = 0 and  y ■ Oslng  these  guidelines,  locations  of  vortices  In  the 

near  wake  can  be  estimated  and  they  are  shown  In  Figure  15.  In  this 
diagram  the  circles  signify  vortex  core  sizes,  and  the  arrows  show  the 
direction  of  circulation.  The  diagram  shows  the  size  of  the  first  vortex 


not  fully  grown  at  inch,  which  could  not  be  established  using  the 

previously  described  criteria  so  the  size  shown  In  the  figure  Is  the 

largest  that  can  be  accommodated  by  the  geometry.  The  dotted  line  Is  the 

— 1/2 

locus  of  the  secondary  maxima  of  u Interpreted  as  the 

boundary  of  the  region  of  secondary  flow  at  the  trailing  edge. 

Also  shown  In  Figure  17  Is  the  streamwlse  distribution  of  the  root- 
mean-square  fluctuating  pressure,  which  has  been  normalized  on  the 
dynamic  head  of  the  free  stream.  It  Is  shown  that  the  wake-induced 
pressure  diminishes  to  one-third  Its  value  at  the  tip  within  only  one 
tralllng-edge  thickness.  These  levels  were  obtained  from  dimensionless 
spectra  such  as  those  In  Figure  15.  In  the  latter  figure  the  frequency 
of  shedding  has  been  made  dimensionless  on  the  cross-wake  distance  be- 
tween filtered  turbulence  maxima  y * y^,  evaluated  at  r ^ ^nd  a 

•'o  •'f  X f 

speed  U , characteristic  of  the  shear- layer  at  the  point  of  separation, 

s g 

This  speed  has  been  defined  by  Roshko  through  the  theory  of  potential 
flow  to  be  related  to  the  static  base-pressure  coefficient  by 


U^/1  - C 

00  p 


The  dimensionless  f lltered-pressure  levels  show  a maximum  at  a dimensionless 

frequency  w y,/U  * 1.  The  bandwidth  of  the  pressure  levels  appears 
S X s 

greater  for  a speed  of  50  than  for  100  ft/sec;  however,  at  both  speeds 

the  bandwldths  of  the  measured  pressures  are  controlled  by  the  filter 

for  which  the  dimensionless  bandwldths  are  Au)y^/U^  « 0.09  and  0.18  for 

100  and  50  ft/sec,  respectively.  Also  apparent  are  the  lower  magnitude 

harmonics  of  the  fundamental  pressure  fluctuations.  For  toy-/U  < 3,  the 

X s 

broadband-pressure  levels  at  the  tip  are  a factor  of  6 to  10  In  excess 
of  those  measured  upstream  at  position  6,  r^  * -13  Inches.  At  the  latter 
position,  spectrum  levels  are  typical  of  those  measured  In  turbulent 
boundary  layers  with  vanishing  streamline  pressure  gradient. 

Correlations  of  the  periodic  pressures  along  the  span  were  determined 
at  a speed  of  ••  100  ft/sec  for  microphone  separations  to  r^  = 4 h. 


m 
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These  correlations  have  been  normalized  on  the  mean-square  pressure,  which 
was  uniform  along  that  portion  of  the  span  used  for  the  measurement. 

This  correlation  is  written 


<p(z,t)  p(z+r  , t+T)> 

R (r  .T) 5 

PP  z 2 ^ 

[p_  (z.t)  p (t,z+r  )] 
z z 


1/2 


(13) 


and  the  measurements  behave  as  A cos  + ct) , where  A and  a are 

separation-dependent  amplitude  and  phase  functions  approaching  unity  and 

zero,  respectively,  as  r^  -+■  0.  The  space-time  correlations  indicated 

that  A was  a slowly  varying  function  of  time  delay  l.  The  correlations 

at  T = 0 produced  the  lateral  coherence  function  A cos  a;  at  time  delays 

T ■ a/oi  , they  produced  the  maximum  correlation  amplitude  A.  The  angle 
m s 

a is  related  to  the  yaw  angle  4>  made  by  the  vortex  filaments  with  the 

trailing  edge  by  the  relationship  a > (o)  r tan  (j>)/U  . The  data  shown  in 

Figure  18  are  for  both  of  these  correlation  functions.  The  scatter  for 

R (r  , T * 0)  was  caused  by  the  slow  variation  of  this  function  with 
PP  z 

time,  reflecting  that  (j>  was  not  constant.  The  phase  of  the  correlation 
a,  in  fact,  varied  between  zero  and  the  absolute  value  reflected  by  the 
data  given  in  Figure  18.  This  was  determined  by  measuring  the  correlation 
with  short  averaging  times.  The  data  indicate  that  the  vortex  filaments 
were  shed  at  yaw  angles  to  (|)  * 17  degrees.  Since  |r  (r  ,0)1  < R (r  ,T  ), 
the  data  of  Figure  18  show  that  for  r^  > 4 h the  lateral  correlation  will 
have  some  negative  values  with  an  absolute  value  less  than  0.2.  A 
correlation  length  defined  as 


i 

z 
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R^^(r  .0)  dr^ 

PP  z z 


(14) 


is  Indicated  by  the  data  to  be  on  the  order  of  = 4 h.  Similarly 
defined  correlation  lengths  have  been  found  for  vortex  streets  behind 
cylinders  to  be  as  low  as  3 and  as  large  as  10  cylinder  diameters,  and 
they  depend  on  Reynolds  number. 

4.2  WAKE  BEHIND  45-DEGREE— ROUNDED- 
BEVELED  TRAILING  EDGE 

Our  discussion  of  this  edge  will  parallel  that  of  the  blunt  edge. 
When  the  trailing  edge  is  unsymmetrically  rounded,  using  the  45-deg  angle 
and  a 5-inch  fairing  radius  of  curvature  as  shown  in  Figure  1,  the 
separation  point  for  the  boundary  layer  on  the  curved  side  becomes  less 
distinct.  This  degree  of  disorder  as  well  as  the  streamwlse  displace- 
ment of  the  free-shear  layer  formation  points  detunes  interaction. 
Characteristics  of  the  boundary  layer  and  near-wake  are  shown  in 
Figure  19.  The  figure  shows  a favorable  pressure  gradient  upstream 
from  tap  11,  followed  by  an  adverse  gradient  downstream  from  this  point 
and  constant  pressure  downstream  from  tap  13.  Observation  of  oil-streak 
patterns  disclosed  concentrations  of  oil  along  spanwise-oriented  lines 
at  taps  13  and  17.  The  former  is  interpreted  as  an  upper  separation 
point;  the  latter,  as  the  stagnation  point  for  the  trailing  edge.  The 
constant  static  pressure  and  the  oil-film  patterns  indicate  that  flow  is 
separated  downstream  from  tap  13. 

The  center  diagram  of  Figure  19  shows  the  mean  velocity  profiles  in 
the  tralling-edge  region.  The  profile  at  position  G is  given  in 
dimensionless  form  in  Appendix  B.  Speeds  measured  at  various  points 
above  the  surface  have  been  made  dimensionless  by  using  a local  approxi- 
mation of  the  free-stream  speed  U^,  which  was  determined  at  each  measure- 
ment station.  The  bulge  in  the  profile  at  tap  11,  showing  a somewhat 
higher  speed  near  the  wall  than  in  the  outer  flow  reflects  the  negative 
pressure  coefficient  at  this  point.  It  is  apparent  that  Equation  (1) 
is  not  satisfied  by  the  data  because  unity  minus  the  square  of  the 
measured  U /U„,  where  U is  the  local  maximum  speed  just  outside  the 
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boundary  layer,  does  not  equal  the  minimum  pressure  coefficient  at  the 
corresponding  point.  The  discrepancy  Is  probably  due  to  an  error  In 
measuring  U^,  determined  at  this  chordwlse  location  about  1.5  Inches  above 
the  surface. 

The  velocity  profiles  show  the  flow  reversal,  which  is  characteristic 
of  regions  of  separated  flow.  Profiles  of  broadband,  rort-mean-square 
velocity  fluctuations  are  also  shown,  and  they  display  local  maxima  at 
distances  from  the  wall  corresponding  to  maxima  In  the  vertical  gradient 
of  mean  velocity  dU/dy.  The  velocity  spectral  densities,  measured  at 
= 100  ft/sec,  disclosed  a periodic  contribution  at  to  ■ to^,  which  is 
Interpreted  as  due  to  a vortex  street.  Profiles  of  narrowband  fluctuating 

velocity,  filtered  In  1/10-octave  bands  at  the  shedding  frequency 

~2  1/2 

0)  ■ 0)^,  u shown  at  the  bottom  of  Figure  19.  Above  the 

surface  of  the  edge  and  In  the  wake  at  downstream  distances  r^  > 1 Inch, 
there  are  maxima  corresponding  to  maxima  In  dU/dy.  Very  near  the  edge, 
there  are  secondary  maxima  which  are  proposed  to  be  associated  with  a 
secondary  vortex  system  similar  to  that  already  proposed  behind  the  blunt 
edge.  In  contrast  to  the  wake  of  the  blunt  edge,  however,  there  Is  an 
asymmetry  In  the  wake  which  persists  at  least  5 Inches  downstream  from 
the  trailing  edge.  In  this  case,  local  maxima  In  filtered  Intensities 
occur  at  the  upper  and  lower  cross  wake  positions  y “ and  y * y^^, 
respectively.  The  characteristic  cross  wake  dimension  is,  then,  y^ 

= ^u  - yr 

A characteristic  shedding  velocity  will  again  be  defined  in 
terms  of  the  pressure  coefficient  -C^  beneath  the  separated  boundary 

® 8 

layer  downstream  from  tap  13.  Following  the  Roshko  definition  we  write 


vA  - C 


(12) 


This  Is  the  same  as  the  separation  velocity  of  the  blunt  edge.  Spectral 


densities  of  the  wake-velocity  fluctuations  are  nondlmenslonallzed 


on  the  wake  variables  y and  in  Figure  20  (the  spacing  y^  will  be 


calculated  later).  The  mean-square  velocity  Is 
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u * I 4>  (to)  d u)  (15 

I uu 
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At  a distance  “ 1.25  inch  downstream  from  the  edge,  and  2y/h  = 1.10 
above  the  tip,  the  dimensionless  spectra  show  the  periodic  contribution 
at  (0  Yf/U  * 1 for  U = 100  and  50  ft/sec.  The  Intensity  of  the 
contribution  Is  substantially  higher  at  100  than  at  50  feet  per  second. 
Velocity  fluctuations  at  ■ 50  ft/sec  and  2y/h  = 0 are  entirely  broad- 
band. A similar  nondlmenslonallzatlon  of  fluctuating  surface  pressures 
on  the  edge  also  shows  a periodic  contribution  at  o)  y^/U  “ 1. 

Figure  21  shows  the  dimensionless  spectra  of  narrow-frequency-band 
(12.5-Hz)  levels  of  mean-square  spectra  at  positions  6 and  G.  For 
“ 100  ft/sec,  the  periodic  contribution  is  greatest  at  position  G 
and  substantially  diminished  at  position  6.  At  50  ft/sec,  the  spectrum 
at  position  G shows  a broadly  peaked  contribution  which  contrasts  the 
weak  peridocity  Indicated  by  the  velocity  spectrum  of  Figure  20.  It 
appears  that  the  pressure  fluctuations  extant  on  the  edge  at  50  ft/sec 
are  determined  more  by  local  separation  than  by  the  vorticity  in  the  wake 
field.  At  100  ft/sec,  the  pressures  generated  by  the  wake  field  are 
dominant . 

Further  corroboration  for  this  hypothesis  is  given  by  the  space- 

time  cross-correlations  between  the  unfiltered  trailing-edge  pressure  at 

position  G and  the  wake-velocity  fluctuations  filtered  in  1/10-octave 

bands  at  f = f^,  similar  to  those  of  the  previous  section.  For 

= 50  ft/sec,  correlations  between  the  pressure  fluctuations  at 

position  G and  the  wake  velocities  along  y = y and  y = y,,  are  shown  in 

li  x» 

Figure  22.  Normalized  on  the  root-mean-square  pressure  at  G and  the 
mean  velocity  U^,  the  spatia.’  correlation  along  the  locus  shows  a 

minimum  of  R (0,r  ,y  ) = -0.04  over  the  edge  whose  absolute  value  is 
never  exceeded  for  r^  > 0.  This  shows  that  instantaneously  the  pressure 
at  position  G is  dominated  locally  by  separation.  The  correlations  also 


show  that  the  shear  layers  along  and  y^  interact  to  form  a coherent 
wake;  however,  the  level  of  coherence  with  the  pressure  at  position  G 
decreases  as  r^  increases.  The  spatial  correlation  indicates  that  the 
formation  of  vortices  from  the  upper  and  lower  shear  layers,  although 
facilitated  by  mutual  interaction  of  the  layers,  does  not  Involve  total 
entrainment  of  one  layer  by  the  other.  Figure  22a  suggests  that  the 
upper  layer  forms  a vortex  first  near  r^  * 1 while  the  lower  layer  forms 
a vortex  near  r =1.5  inch.  This  indicates  a lag  on  the  order  of  one- 
eighth  wavelength.  Further  downstream,  the  mutual  interaction  brings 
the  vortices  into  a more  coherent  alignment. 

This  interpretation  of  the  spatial  correlation  is  also  supported 
by  the  growth  of  the  root-mean-square — f iltered-velocity  fluctuations 
of  Figure  19.  The  local  maximum  values  of  filtered  disturbances  for  both 
shear  layers  are  shown  as  a function  of  r^  in  Figure  22b.  Along  y “ 
the  disturbances  filtered  at  u>  = grow  to  a maximum  at  a distance 
between  1.5  inch  < r <3  inch,  downstream  but  along  y = yA  the 

X JC 

disturbances  grow  to  a maximum  at  r^  = 1 inch.  The  growth  and  streamwlse 
coherence  of  the  disturbances  shown  by  Figures  22a  and  22b  indicate  that 
the  development  of  regularly  spaced  vortices  is  weakened  by  the  relatively 
Independent  behavior  of  the  upper  and  lower  shear  layers  as  compared  to 
that  behavior  of  those  shed  from  the  blunt  edge. 

The  growth  and  decay  of  a given  eddy  as  it  moves  downstream  from 

position  G as  indicated  by  the  moving  eddy  correlation  is  shown  in 

Figures  22c  and  22d.  The  time  T has  been  defined  from  the  maximum 

m 

magnitude  of  the  space-time  correlation  as 


|r  (r  ,T  )|  > |R  (r  ,T)| 
' pu  X m ' ' pu  X ' 


where  r^  is  fixed.  Along  the  trajectory  of  the  upper  shear  layer 
y = y^,  the  absolute  value  of  the  correlation  grows  to  give  a minimum 
at  r^  = 1.5  Inch,  Rp^(1.5  in.,  T^)  = -0.052.  Along  the  lower  trajectory. 


the  value  of  the  space-time  correlation  is  taken  at  the  time  delay  of 
maximum  positive  correlation  since  we  are  interested  in  a contribution 
180  deg  out  of  phase  with  the  observed  eddy  at  y ■ y^.  This  phase 
reversal  is  dictated  by  the  necessity  for  the  vortices  of  each  shear 
layer  to  be  of  opposite  sign.  The  correlations  in  the  upper  and  lower 
shear  layers  do  not  become  matched  until  r^  > 2 Inches;  the  time  delays 
become  nearly  equal  at  this  distance. 

The  behavior  of  the  near  wake  for  the  45-deg  round  edge  is  strongly 
dependent  on  Reynolds  number  as  can  be  seen  now  by  reference  to  1 guie  23. 
The  correlation  measurements  were  repeated  for  = 100  feet  per  second. 

In  contrast  to  the  case  of  low  speed,  the  disturbances  from  both  layers 
are  seen  to  grow  to  maxima  at  1.5  inch  < r^  < 2 Inches  with  the  spatial 
correlation  increasing  dramatically  with  distance  downstream  from  the 
tip  of  the  trailing  edge.  In  this  case  < 0.1,  and  r^  > 0, 

compared  to  R (r  ,0)  < 0.023  for  r >0.  The  contribution  to  the  total 

pu  X - X 

pressure  fluctuation  by  the  local  separation  is  notably  reduced,  compared 
to  the  contribution  at  U =50  ft/sec  as  shown  by  the  lower  correlation 
magnitude  f>)r  r^  < 0.  It  is  worthwhile  to  note  that  the  root-mean- 

I square  pressure  that  was  used  to  normalize  the  correlation  is  the  broadband 

~2  1/2 

pressure,  p . Reference  to  Figure  21  shows  that  while  the  periodic 

! pressures  are  higher  and  dominate  the  signal  at  = 100  ft/sec,  the 

broadband  pressures  are  dominant  at  = 50  feet  per  second.  For  this 

~2  1/2  _2  ~2  1/2 
’ reason,  p /q^  = 1.13  x 10  at  100  ft/sec;  however,  p /q^ 

-2 

I = 2.6  X 10  at  50  feet  per  second.  Thus  it  appears  that  even  though  a 

development  of  the  periodic  wake  is  weak  at  low  speeds,  the  continuous 

j!  spectrum  pressures  due  to  the  local  separation  are  Important.  The 

j moving-eddy  correlation  ^oo  “ 1®®  ft/sec  not  only  reflects 

the  pronounced  wake  development  but  the  values  of  show  small 

At  /Ar  for  small  r . This  is  behavior  quite  similar  to  that  observed 
J m X X 

for  the  blunt  edge.  Also,  the  values  of  R (r  ,T  ) for  the  upper  and 

pu  X m 

lower  shear  layers  are  more  nearly  matched  at  = 100  feet  per  second. 
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A wake  structure  for  U = 100  ft/sec  as  deduced  from  the  measure- 

00 

merits  is  shown  in  Figure  24.  As  in  Section  4.1,  cores  of  the  vortices 
are  assumed  to  be  bounded  by  the  lines 


y > y > 
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with  a streamwise  spacing  given  by  R Figure  23.  The  dotted 

line  near  the  trailing  edge  circumscribes  both  the  region  of  flow  reversal 

1/2 

as  well  as  the  locus  of  secondary  maxima  in  u (f^)  as  shown  in 

Figure  19  for  small  r . Root-mean-square  filtered  pressures 

-y  1/2  ^ 

P (f  )/q„  are  shown  to  have  a maximum  at  position  G(r  ■=  -1.1). 

At  the  downstream  position  (r^  = -0.4)  the  periodic  contributions  are 
missing;  this  position  is  nearly  coincident  with  the  stagnation  point 
r ■ -0.6.  The  convection  velocities 

X 
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as  determined  from  Figure  23  show  a maximum  at  r^  = 0.8  along  y * y^. 

Here  the  velocity  disturbance  is  approximately  180  deg  out  of  phase  with 
the  fluctuating  pressure  at  G,  and  its  phase  slowly  varies  as  the  eddy 
moves  downstream.  Along  the  lower  layer  y = an  opposite  behavior 
appears;  the  disturbance  convection  velocity  decreases  slightly.  This 
contrasting  behavior  could  be  caused  by  an  uncertainty  in  determining 
the  peak  in  the  space-time  correlation  and  would  be  discounted  altogether, 
except  that  Hanson^^  made  a similar  observation  in  the  wake  behind  a 
notched,  lifting  airfoil.  He  observed  convection  velocities  on  the 
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order  of  1.7  U^.  The  current  measur&nent  Indicates  that  the  upper 
disturbance  speeds  up  briefly,  while  the  lower  disturbance  Is  convected 
at  a lower  speed.  The  formation  length  defined  as  the  distance  from  the 
stagnation  point  on  the  foil  to  the  value  of  r^  for  the  maximum  in  the 
filtered  velocity  fluctuation,  is  taken  to  be  = 2.3  Inches.  This 
nearly  corresponds  to  the  location  of  the  minimum  shear  layer  spacing 

" ^u  - yji  “ yf 

4.3  CHORDWISE  DISTRIBUTION  OF 
FLUCTUATING  PRESSURE 

Figures  17  and  24  show  that  the  root-mean-square  pressures  on  each 
edge  decrease  with  Increasing  upstream  distance  from  the  extremity  of 
the  edge.  In  Figure  25,  mean-square  filtered,  12.5-Hz  pressures*  for 
each  edge  are  shown  together  as  a function  of  dls'.ance  from  the  stag- 
nation point  (x-x  ).  The  measurements  were  made  on  each  edge  for  two 
8 

speeds,  and  the  pressures,  normalized  on  the  dynamic  head,  have  been 
grouped  for  a given  edge  shape.  The  chordwlse  variation  for  the  blunt 
edge  is  less  dependent  on  speed  than  for  the  45-deg  round  edge.  This 
difference  could  be  due  to  the  occurrence  of  three  dimensionality 
associated  with  a finite,  spanwise  correlation  of  the  vortex  streets. 

In  this  regard,  recall  that  Figure  18  gives  a spanwise  correlation  length 
« 4 h ■ 8 inch.  The  data  of  Figure  25  show  that  forward  of  the  blunt 
edge  for  (x-x^)  > 10  inches,  the  rate  of  decrease  of  pressure  with 
further  distance  is  more  pronounced.  This  suggests  a three  dimensionality 
effect  on  the  chordwlse  pressure  distribution  which  may  be  speed 
dependent.  On  the  45-deg  round  edge,  the  less  rapid  reduction  in 
fluctuating  pressure  with  distance  at  U^^  ■■  200  ft/sec  could  be  attributed 
to  greater  spanwise  coherence  at  the  greater  speed.  Unfortunately, 
spanwise  correlation  measurements  were  not  made  for  this  edge. 

A dimensional  argument  for  scaling  these  pressure  distributions  is 
derived  from  the  estimated  wake  structures  in  Figures  17  and  24.  For 


Filter  bandwidth  was  greater  than  bandwidths  of  pressures. 
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a vortex  situated  at  downstream  from  the  stagnation  point,  the  local 
pressure  exerted  on  the  strut  at  x Is  assumed  to  be  determined  In  a 
general  functional  form  by 


P(x,x  ) 
s 


(=?) 


where  U Is  the  shedding  velocity  defined  earlier  on  the  basis  of  the 
static  pressure  coefficient  in  the  trailing-edge  separation 
region 

x-Xg  Is  the  position  of  the  vortex  relative  to  the  edge 

K Is  the  strength  of  the  vortex. 

The  Influence  function  g((x-Xg)/H^)  specifies  that  the  pressure  at  a 
point  Is  governed  only  by  the  distance  to  the  vortex.  The  function  is, 

In  general,  dependent  on  the  shape  of  the  edge.  For  a periodic  vortex 
sheet,  the  Individual  vortices  are  correlated  with  a specific  spatial 
relationship  so  that  the  formation  length,  which  locates  the  first  vortex 
pair  relative  to  the  edge,  is  a dominant  space  variable. 

The  characteristic  strengths  of  the  vortices  In  the  wake  are 

determined  by  the  core  radii  r given  in  Figures  17  and  24  and  the  root- 

° -2  1/2 

mean-square  filtered  Intensities  u (f  ) of  Figures  16  and  23.  The 

s 

second  vortex  is  considered  to  be  characteristic  of  the  wake  because  for 

each  edge  It  is  the  first  vortex  immediately  downstream  from  formation. 

— 1/2 

In  the  case  of  the  blunt  edge,  ^r^/h  = 0.8,  and  u*"  (f^)  •*  0.29 

so  that  we  approximate 

-j  1/2 

K=/2»27ru  (f)*r 

so  / -I  o \ 


K = 0.52  U h 
o 


as  the  typical  vortex  strength.  Greenway  and  Wood^  have  estimated  that 


0.5  < </U  h < 1.2 


from  their  flow  visualization  of  blunt-edge  flows  at  chordwlse  Reynolds 


numbers  of  1.4  x 10  . At  the  decade  higher  Reynolds  number  of  the 


current  measurements,  taking  U = U , the  vortex  strength  appears  to  be 
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slightly  reduced.  For  the  45-deg  round  edge,  4r  /h  = 0.5,  and  u^  (f  ) 

o s 


0.11  U^,  so  that  K = 0.12  U^h.  Thus  rounding  the  trailing  edge  not  only 


decorrelated  the  wake  but  also  reduced  the  coherent  vortex  strength. 

The  normalization  of  Equation  (17)  describes  a universal  distribution 
function.  To  test  its  validity,  the  data  of  Figure  25  have  been 


replotted  in  a dimensionless  form,  utilizing  as  a length  scale  in 


Figure  26.  The  pressures  appear  to  be  well  defined  by  these  variables — 
some  points  for  200  ft/sec  have  been  Ignored  in  this  representation. 


Using  reciprocals  of  integral  powers  of  distance  we  find  that  a series 


-1  -3 

combination  of  (x-x  ) and  (x-x  ) would  probably  describe  the 

s s 


dimensionless  distribution  over  a significant  range  of  the  variable 


(x-Xg)/5.j. 


4.4  VORTEX  SHEDDING  FREQUENCY; 
UNIVERSAL  STROUHAL  NUMBER 


The  vortex  shedding  frequency  has  been  shown  by  Figures  15  and  20 


or  21  to  scale  on  the  shedding  speed  U^  and  the  near-wake  length  scale 


y, . An  alternative  nondlmenslonallzatlon  which  has  been  suggested  by 
^ 10 

Hanson  Involves  the  momentum  thickness  of  the  wake  at  formation 


0^  where 
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His  data  for  lifting  airfoils  show  that  the  Strouhal  number 


f e,/u 

s f s 


is  a universal  function  of  the  Reynolds  number  9^/v.  However,  the 
data  for  the  current  edges  only  partially  support  this  definition  of 
Strouhal  number. 

Figure  27  shows  the  Hanson  results  as  dots  together  with  the  points 
which  apply  to  the  current  experiments.  The  curve  drawn  through  the 
points  follows  the  empirical  equation 


Sq  = 0.0728  [Nj^  - 1038]  (IS 

which  was  determined  by  Hanson  and  which  is  similar  to  that  determined 

g 

by  Roshko.  The  dimensionless  frequency  for  the  shedding  from  the  blunt 
trailing  edge  is  only  one-half  that  given  by  the  45-deg,  round  edge  and 
by  Equation  (19).  Apparently,  the  momentum  thickness  of  the  wake  at 
formation  is  not  a universally  applicable  length  scale  for  the  shedding 
frequency.  A distinguishing  feature  of  blunt  edge,  compared  to  the 
45-deg,  round  edge  and  the  edges  of  the  Hanson  study,  seems  to  be  the 
base-pressure  coefficient  which  is  (-C  )_  * 0.51  in  the  first  case  and 
less  than  (-C  )„  = 0.3  in  all  of  the  later  cases.  The  Bearman"^  and 

10  P ® 

Hanson  investigations  with  splitter  plates  have  shown  that  an  increase 
in  base-pressure  coefficient  accompanies  a reduction  in  formation 
length.  Also,  a relatively  low  base-pressure  coefficient  represents 
more  complete  pressure  recovery  by  the  flow  and  this  in  turn  implies 
that  less  momentum  has  been  transferred  to  the  wake  in  the  form  of 
fluid  unsteadiness  at  the  trailing  edge.  The  loss  of  momentum  from  the 
mean  flow  downstream  from  the  edge,  which  determines  6^,  must  be  balanced 
by  the  Increase  of  unsteady  fluid  momentum.  The  growth  of  the  regular 
disturbances  associated  with  periodic  vortex  streets  is  undoubtedly 
accelerated  by  the  presence  of  this  unsteady  motion;  however,  it  is 
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probably  not  fully  governed  by  it.  Thus  the  wake  momentum  thickness, 
while  being  representative  of  the  overall  loss  of  momentum  from  the 
mean  flow,  may  not  be  expressive  of  the  dynamical  conditions  necessary 
for  the  formation  of  a periodic  vortex  street.  The  conditions  favorable 
for  the  formation  of  regular  disturbances  might  be  expected  to  be  governed 
by  6^  rather  than  by  a geometric  dimension  of  the  trailing  edge,  e.g., 
the  lateral  distance  between  the  flow  separation  points,  if  the  vortices 
are  formed  far  enough  downstream  from  the  trailing  edge  that  extraneous 
flow  disturbances  associated  with  the  immediate  separation  region  decay. 
This  is  probably  not  the  case  for  very  blunt  trailing  edges  with  high 
base-pressure  coefficients  where  vortices  are  formed  close  to  the  edge 
but  may  be  possible  for  smoother  edges  with  lower  C 

^b 

The  formation  of  the  regular  vortex  street  from  the  shear  layers 
leaving  the  shedding  body  will  depend  on  the  level  of  vortlclty  associated 
with  each  layer  as  it  is  shed  as  well  as  on  the  cross-wake  distance  be- 
tween the  layers.  The  following  discussion  will  show  the  motivation  for 
suggesting  y^  as  an  alternative  to  6^  as  a length  scale.  Some  theoretical 

work  has  been  performed  to  examine  the  influences  of  these  factors. 

22 

Abernathy  and  Kronauer  (1962)  have  performed  calculations  of  the 

growth  of  disturbances  in  idealized  parallel  shear  layers.  The  degree 

of  shear-layer  interaction  for  a given  streamwlse  wavelength  of  disturbance 

was  shown  to  depend  on  the  separation  of  the  layers.  A breakdown  in  the 

shear  layers,  due  to  instability  to  small  'isturbances,  was  shown  to 

occur  with  the  resultant  formation  of  idenf ' 'illy  structured  concentrations 

of  vortlclty.  The  spatial  structure  of  these  disturbances,  which 

resembled  that  of  a von  Karman  vortex  street,  motivated  Abernathy  and 
22 

Kronauer  to  speculate  that  there  is  a unique  relationship  between  the 
wavelength  of  the  initial  disturbance  and  shear-layer  separation 
distance  which  is  most  likely  to  generate  a periodic  vortex  street. 
Physically,  thin  sheets  of  vortlclty,  such  as  those  which  are  assumed 


Abernathy,  F.H.  and  R.E.  Kronauer,  "The  Formation  of  Vortex  Streets," 
Journal  of  Fluid  Mechanics,  Vol.  13,  pp.  1-20  (1962). 


in  the  idealized  calculations  of  Abernathy  and  Kronauer,  appear  to  have 

existed  only  in  the  immediate  vicinity  of  the  blunt  trailing  edge  as 

shown  in  Figure  14.  They  are  not  apparent  at  all  behind  the  45-deg  round 

edge  (Figure  19)  where  the  shear  layers  are  relatively  thick.  To 

examine  the  stability  of  disturbances  in  wakes  which  are  characterized 

by  thick  shear  layers  similar  to  those  of  the  round  edge  (Figure  19)  a 

23 

stability  analysis  has  been  performed  by  Sato  and  Kurlkl  (1961)  . They 

found  that  at  large  Reynolds  numbers,  small  disturbances  of  a wide  range 

of  wavelengths  could  be  amplified  and  that  the  rate  of  the  amplification 

was  determined  by  the  inflection  points  in  the  mean  velocity  profile. 

The  calculations  were  verified  by  measuring  disturbances  in  the  wake  of 

a thin  flat  plate.  The  distance  between  the  Inflection  points  of  the 
23 

Sato  and  Kurlkl  profile  was  essentially  analogous  to  the  shear  layer 

22 

separation  of  Abernathy  and  Kronauer.  Also,  the  amplitudes  of  the 

23 

g-  owing  disturbances  were  observed  by  the  Sato  and  Kuriki  combined 
measurement  and  calculation  to  be  maximum  at  the  inflection  points  of 
the  mean  velocity  profile. 

The  wake  stability  calculations  described  in  Appendix  C were 
performed  for  a variety  of  profile  shapes  of  the  general  form  of  those 
in  Figure  19.  They  have  further  demonstrated  that  the  frequency  of  the 
least  stable  disturbances  in  each  case  depends  more  on  separation  of  the 
Inflection  points  than  on  the  momentum  thickness  of  the  wake.  The  momentum 
thickness  of  each  profile  is  shown  to  be  not  uniquely  related  to  the 
shear-layer  separation  for  the  set  of  profiles  considered. 

These  considerations  motivate  the  selection  of  y^  as  an  alternative 
to  0^  as  a near  wake-length  scale  for  the  spectra  of  Figures  15  and  20 
or  21.  Use  of  as  a speed  scale  is  motivated  by  the  view  that  it  is 
representative  of  the  speed  of  the  shear  layer  as  it  leaves  the  shedding 
body.  Further  evidence  from  past  investigations  for  the  use  of  this 
frequency  scaling  is  scant  because  y^  is  not  available  In  the  literature. 

23 

Sato,  H.  and  K.  Kurlkl,  "The  Mechanism  of  Transition  in  the  Wake  of  a 
Thin  Flat  Plate  Placed  Parallel  to  a Uniform  Flow,"  Journal  of  Fluid 
Mechanics,  Vol.  11,  p.  321  (1961). 
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Neither  the  formation  length  nor  were  determined  in  earlier  investi- 
gations other  than  those  of  Bearman^  and  of  Hanson. The  measurements 
of  Hanson^^  which  are  among  the  most  extensive  in  this  regard  are  not 
complete  enough  to  determine  very  precisely.  Table  3 shows  the  best 
estimates  of  U)  y./U  which  can  be  extracted  from  data  by  Bearman  and 
Hanson. 


TABLE  3 - STROUHAL  NUMBERS 


Edge  Shape 

0)  y-/U 
s f s 

Position  of  y^ 

Blunt^ 

0.97 

1.0 

Blunt  with  splitter  plate, ^ i/h=l 

0.92 

1.0 

Notched  lifting^®  (N^^  = 6900) 

®0 

0.78 

0.3 

Notched  llftlng^^  (Nj^  * 4800) 

0.83 

0.5  Si^ 

Splitter, )l/h=l  (Nj^  » 5930) 

®0 

0.91 

0.4 

Splitter, Jl/h-1  (Nj^  - 9800) 

®0 

0.94 

0.5 

The  wake  position  downstream  from  the  trailing  edge  at  which  y^  is  esti- 
mated has  also  been  listed.  I,t  is  pertinent  to  note  that  the  lowest 

value  of  u)  y,/U  occurs  for  the  case  of  y,  measured  only  0.3  down- 
s f s ■'f  t 

stream  from  the  edge.  Because  of  wake  spreading  we  would  expect  the 
actual  value  of  y^  to  be  somewhat  greater  than  that  used  here.  This 
nondlmenslonallzatlon  appears  to  scale  the  shedding  frequency  within 
20  percent. 

Bearman^  has  suggested  an  alternative  scaling  based  on  far-wake 

9 

variables.  This  scaling  is  based  on  the  Kronauer  stability  criterion 
which  states  that  the  stable  cross-wake  spacing  b of  the  vortices  in  a 
wake  with  a periodic  vortex  street  is  of  such  a nature  that  a minimum 
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drag  is  exerted  on  the  body.  The  Bearman  universal  Strouhal  number 

S„  = f b/U  was  determined  from  the  far-wake  data  of  both  the  blunt 
B s s 

and  45-deg-round  edges  using  his  calculation  procedure.  The  values  are 
shown  in  Figure  28  as  a function  of  shedding  velocity  U^/U^;  agreement 
with  the  trend  of  earlier  data  continues  to  support  the  Bearman  definition. 


4.5  SUMMARY:  PERIODIC  VORTEX  STREETS 

The  temporally  periodic  pressures  that  are  generated  on  a strut  by 
the  formation  of  a von  Karman  vortex  street  have  been  quantitatively 
related  to  the  position  of  the  first  fully  formed  vortex  in  the  wake  as 
well  as  to  the  strengths  of  the  vortices.  This  relationship  gives 
the  chordwise  variation  of  the  magnitude  of  the  pressure  as  a function 
of  the  distance  upstream  from  the  edge  stagnation  point  expressed  in 
multiples  of  the  vortex  formation  length.  This  similarity  function 
appears  to  be  independent  of  both  the  speed  and  the  edge  geometry  over 
the  range  of  variables  considered.  Some  influence  of  the  spanwise 
coherence  of  the  vortex  street  on  the  chordwise  behavior  is  indicated. 

A spanwise  coherence  length  of  approximately  four  edge  thicknesses  was 
observed  for  the  blunt  edge.  For  a fixed  formation  length,  the  magni- 
tude of  the  pressure  is  directly  proportional  to  the  strengths  of  the 
vortices.  Development  of  the  vortex  street  was  suppressed  at  low 
speeds  behind  the  45-deg-round  edge;  disturbances  decayed  as  they  were 
convected  downstream  from  the  edge.  In  this  case  periodic  pressures 
could  not  be  detected  on  the  strut.  At  higher  speeds  the  spatial 
wavelike  growth  of  disturbances,  typical  of  the  blunt-edge  wake,  was 
detected,  and  periodic  pressures  were  generated.  At  the  lower  speeds, 
however,  even  though  the  periodic  disturbances  were  suppressed,  high- 
level,  continuous-spectrum  pressures  were  generated  by  the  local 
boundary-layer  separation. 

The  frequency  of  vortex  shedding  was  found  to  be  given  by  two 
Strouhal  numbers  which  are  based  on  either  the  near-  or  far-wake 
variables.  The  near-wake  momentum  thickness  was  verified  as  a 
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characteristic  length  for  the  shapes  with  low  base-pressure  coefficients. 
The  estimated  shear  layer  separation  at  the  end  of  the  formation  region 

was  found  to  be  a characteristic  length  with  more  general  application. 

7 9 

The  Bearman  far-wake  Strouhal  number,  based  on  the  Kronauer  stability 

criterion,  was  found  to  agree  with  the  current  data. 

5.  SEPARATED  FLOW  AT  BEVELED  TRAILING  EDGES 

The  remainder  of  the  report  will  be  concerned  with  flow  near  the 
nonsinging  edges  of  Figure  1.  In  these  cases  the  boundary  layers 
incident  on  the  edges  experience  adverse  pressure  gradients  before 
being  shed  by  the  strut.  For  nonsymmetrlcal  edges,  the  flow  on  one 
side  of  the  strut  separates  far  enough  upstreata  from  the  apex,  or  tip, 
of  the  trailing  edge  that  periodic  vortex  street  are  not  formed.  This 
section  describes  the  mean  and  unsteady  velocities  that  occur  on  two 
trailing  edges;  Sectloa  6 will  describe  the  characteristics  of 
fluctuating  pressures  that  are  generated  for  each  case  in  the  separation 
process.  The  purpose  of  these  sections  is  to  account  for  the  separated 
flow  field  which  has  been  shown  in  the  last  section  to  Induce  high- 
level,  continuous-spectrum  pressures. 

5.1  TURBULENT  FLOW  NEAR  25-DEGREE— 

KNUCKLE-BEVELED  EDGE 

As  the  boundary  layer  approaches  the  edge,  it  first  experiences  a 
favorable  pressure  gradient,  followed  by  an  adverse  gradient  just  down- 
stream from  the  knuckle.  This  is  shown  in  Figure  29  for  the  two  tunnel 
speeds  U^^j  ■ 50  and  100  feet  per  second.  The  pressure  recovery  appears 
to  be  more  apparent  at  the  higher  speed  since  the  pressure  coefficients 
are  less  negative  at  100  than  at  50  feet  per  second.  Boundary-layer 
characteristics  were  obtained  at  the  positions  A through  F on  the 
surface  of  the  edge  as  well  as  in  the  near  wake  downstream  from  the 
edge.  These  positions  are  shown  for  reference  in  Figure  29.  Before 
discussing  the  measurements  it  is  appropriate  to  mention  that  an  oil- 
streak  experiment — identical  to  that  described  in  Sections  3 and  A — 
showed  separation  at  position  E.  Oil  accumulated  in  a streak 


one-tenth  of  an  inch  in  width  at  the  downstream  edge  of  the  knuckle.  The 
streaks  were  stagnant  downstream  from  position  F.  This  indicates  a 
steady  separation  at  the  knuckle  followed  by  a large-scale  and 
intermittent  separation  zone  further  downstream. 

Results  of  a survey  of  the  boundary  layer  on  the  edge,  performed  at 
> 100  ft/sec,  are  shown  in  Figure  29.  Expanded  and  normalized 
versions  for  speeds  of  50  and  100  ft/sec  are  shown  in  Appendix  B. 

Table  1 gives  thickness  parameters  for  the  boundary  layer.  At 
position  E,  which  is  one-tenth  of  an  inch  downstream  from  the  knuckle, 
no  evidence  of  the  flow  reversal  that  accompanies  separation  could  be 
found.  A single  hot  wire  anemometer  probe  is  expected  to  sense  the 
absolute  value  of  the  mean  velocity  in  the  boundary  layer;  however,  a 

d|u| 

sign  change  in  the  vertical  gradient  , where  |u|  is  the  magnitude 

of  the  mean  velocity  in  the  boundary  layer,  should  indicate  a reversal 
in  flow.  This  sign  change  was  not  detected  at  any  location  on  the 
trailing  edge.  The  profile  at  position  F shows  a vanishing  velocity 
gradient  near  the  wall  and  that  at  position  H shows  a nearly  vanishing 
gradient  at  the  tip  of  the  edge.  The  linearized,  broadband,  root-mean- 
square  turbulence  intensities  of  the  boundary  layers  are  also  shown  as 
dotted  lines  for  each  location.  These  profiles  show  two  important 
features.  The  first  is  the  display  of  local  maxima  in  turbulence  in- 
tensity which  coincide  with  maxima  in  the  vertical  gradients  in  mean 

“2 

velocity.  These  local  maxima  show  absolute  maxima  in  excess  of  u / 

The  intensity  at  position  I is 


= 0.23,  which  occurs  near  position  H. 

comparable  to  that  at  position  H.  The  second  feature  is  that  near  the 

1/2 

M f 

wall;  where  dU/dy  = 0,  the  gradient 


(^) 


also  vanishes.  Further- 


more, at  position  F the  root-mean-square  turbulence  intensity  in  this 
region  near  the  wall  is  approximately  one  half  the  measured  local  mean 
velocity.  Peak  fluctuations  of  velocity,  therefore,  probably  exceed 


AO 


the  mean  velocity  near  the  wall.  Lines  of  constant  speed,  which  describe 
streamlines,  have  been  drawn  In  Figure  29.  The  dashed  line  encloses 
the  region  near  the  wall  for  which  dU/dy  = 0. 

The  short  dashed  lines  In  the  drawings  of  the  mean  velocity  profiles 

at  positions  F and  H reflect  the  interpretation  that  the  flow  near  the 

wall  Is  oscillatory.  The  nature  of  this  oscillation  Is  shown  by  reference 

to  the  space-time  correlation  of  the  streamwlse  component  of  velocity  at 

position  F.  A pair  of  anemometer  probes  was  positioned  a distance  n 

above  and  in  the  plane  with  the  surface  of  the  strut.  In  the  plane  of 

the  edge  surface,  therefore,  the  probes  were  separated  In  the  streamwlse 

direction  a distance  r *>0.5  Inch  and  In  the  lateral  direction  a distance 

s 

r^  “ 0.13  inch.  Since  the  probes  were  at  a common  height  from  the  edge, 
they  were  nominally  separated  a distance  r^  * 0.  Correlations  were 
measured  at  various  distances  above  the  wall  as  the  probes  were  traversed 
together  along  the  y axis.  The  correlations  were  normalized  as 


p ; T'y  - u(x,t)  u(x+r,fK) 

-r  ^ 1/2 

[u  (x,t)u  (x+r,t)] 


(20) 


Coordinates  x pertain  to, the  location  of  the  reference  probe.  The 

separation  coordinates  r are  measured  In  the  plane  of  the  edge.  The 
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interpretation  of  Favre,  Gaviglio,  and  Dumas  (1957)  ’ is  given  to  the 

correlations.  The  results  are  shown  In  Figure  30  tor  reference  probe 
distances  y * 0.1,  0.25,  0.15,  1.0  Inch  from  the  wall  at  position  F and 
for  “ 50  feet  per  second.  Reference  to  Figure  29*  shows  that  the 
position  y - 0.5  nearly  corresponds  to  the  local  maximum  In  turbulence 
intensity  at  position  F.  At  a greater  distance  above  the  edge,  y = 6 
= 1.0  Inch,  the  crossy  correlation  shows  a maximum  at  a time  delay 
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Favre,  A.J.  et  al.,  "Space  Time  Double  Correlations  in  a Turbulent 
Boundary  Layer,"  Journal  of  Fluid  Mechanics,  Vol.  2,  p.  313  (1957). 
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Favre,  A.J.  et  al.,  "Further  Space-Time  Correlations  of  Velocity 
in  a Turbulence  Boundary  Layer,"  Journal  of  Fluid  Mechanics,  Vol.  3,  p.  344 
(1957). 

Velocity  profiles  for  * 50  ft/sec  and  for  and  for  = 100  ft/sec 
are  similar. 
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t * 0.7  X 10  seconds.  This  defines  a convection  velocity 
c 

77^  -/r  /T  - 1.2 
Uoo  \ s c;  » 

That  U > U„  is  of  no  concern,  since  the  direction  of  the  vector  or 

Q 00  ' 

mean  flow  probably  does  not  coincide  with  the  direction  of  the  vector 
of  separation  so  that  the  measured  convection  velocity  deviates  from 
the  actual  convection  velocity,  e.g.,  V'  by 


(21) 


«=  u'/cos  e 


(22) 

t. 


where  6 is  the  angle  between  the  vectors  Just  mentioned.  At  the  maximum 
of  turbulence  intensity  y * 0. 5 inch,  maximum  correlation  is  reduced 
considerably,  and  it  occurs  at  a larger  time  delay.  This  is  consistent 
with  slower  convection  closer  to  the  wall.  Still  closer  to  the  wall  at 
y = 0.25  inch,  corresponding  to  the  outer  limit  of  the  region  where 
dU/dy  * 0,  the  maximum  correlation  is  increased,  and  it  occurs  at  T ■ 0. 

This  indicates  that  the  fluid  disturbances  here  occur  simultaneously  at 
both  probes.  Near  the  wall,  at  y * 0.1  inch,  the  maximum  correlation 
Increases  in  magnitude,  yet  it  remains  at  T * 0.  The  correlations 
near  the  wall  are  marked  by  a skewness  toward  negative  time  delays. 

This  lack  of  symmetry  in  correlation  can  be  explained  by  suggesting  that 
two  velocity  components  exist  in  the  flow  at  this  position.  The 

strongest  component,  with  the  largest  correlation  at  T ■ 0,  is  not  / 

convected;  a weaker  component,  indicated  by  a reduced  maximum  correlation 
value  at  T < 0,  is  convected  upstream  or  opposite  to  the  direction  of 
mean  flow. 

Correlations  between  probes  positioned  vertically,  r = 0 and 

s 

r^  ^ 0,  are  shown  in  Figure  31.  The  reference  probe  is  set  at  y = 0.1  inch 

above  the  edge.  For  r^  = 0.22  inch,  the  space-time  correlation  at 

T = 0 is  positive.  For  r^  * 0.44  inch,  the  correlation  at  T = 0 is 

small  and  negative.  At  this  separation  one  probe  is  situated  0.1  inch  I 


mHBM 


above  the  wall,  and  the  other  is  set  at  the  point  of  maximum  turbulence 
intensity.  The  small  and  negative  correlation  at  x ■ 0 shows  that  a 
large  and  positive — directed  downstream — velocity  disturbance  in  the 
shear  layer  is  weakly  associated  with  a small  and  negative — directed 
upstream — velocity  disturbance  near  the  wall.  The  Increased  negative 
correlation  magnitudes  at  X ■ -12  msec  show  that  this  positive  velocity 
disturbance  in  the  shear  layer  is  followed  by  a negative — directed 
upstream — disturbance  of  moderate  magnitude.  The  velocities  in  the 
outer  flow  and  those  near  the  surface  are  in  quadrature;  a streamline- 
directed  velocity  in  the  outer  flow  is  correlated  with  a deceleration 


of  the  inner  flow. 


Figure  32  shows  spectral  densities  of  the  fluctuating  velocities*  in 
the  region  of  the  trailing  edge,  normalized  on  the  total  mean-square 
velocity  as  described  by  Equation  (3).  At  position  F,  the  dimensionless 
spectra  appear  to  be  similar  at  all  distances  from  the  wall.  This 
suggests  that  eddies  are  convected  at  a velocity  proportional  and  nearly 
equal  to  U(y),  since  a wave  number  U)/U(y)  can  be  used  to  describe  the 
data.  Furthermore,  this  wave  number  can  be  nondlmenslonallzed  on  5^  to 
give  dimensionless  spectra  which  are  typical  of  those  measured  in  .urbu- 
lent  boundary  layers  on  flat  surfaces.  This  can  be  seen  by  comparing 
to  spectra  obtained  at  position  B.  Also,  we  see  that  the  normalization 
applies  for  two  velocities. 

For  distances  above  the  wall  greater  than  0.3  £,  this  normalization 
of  the  spectra  can  be  expected  because  the  space-time  correlations 
(Figure  30)  show  that  the  turbulence  is  convected  at  velocities  which 
are  comparable  to  those  in  the  boundary  layer.  Closer  to  the  wall, 
the  correlations  in  Figure  30  show  that  for  the  most  part,  the  dis- 
turbances at  two  streamwlse-separated  locations  are  recorded 


It  is  to  be  noted  that  these  spectra  as  well  as  those  in  Figure  33  were 
not  all  measured  in  the  Anecholc  Flow  Facility.  Some  were  measured  in 
another,  smaller  facility  during  preliminary  stages  of  this  work.  In  the 
smaller  facility  slight  blockage  effects  caused  thinner  boundary  layers 
than  those  under  discussion;  however,  the  dimensionless  mean-  and  turbulent- 
/eloclty  profiles  were  basically  Identical  in  the  two  facilities.  Thus 
the  dimensionless  spectra  from  those  measurements  are  Included  here.  Where 
spectrum  measurements  in  the  two  facilities  could  be  made,  the  dimensionless 
data  agreed. 


simultaneously  rather  than  being  convected.  The  convection  Indicated 
by  the  existence  of  the  dimensionless  spectra  for  y < 0.3  6 apparently 
plays  only  a small  part  In  the  near-wall  turbulence  dynamics.  At  these 
distances  from  the  wall,  the  turbulence  Is  distributed  over  wave  numbers 
which  are  large  compared  to  a displacement  thickness. 

Also  apparent  In  the  spectra  of  Figure  32  are  certain  peaks  at 
different  low  dimensionless  frequencies  which  actually  occur  at  the  same 
low  frequency.  Normalized  on  the  free-stream  velocity,  the  peaks  occur 
at  o)6*/U^  = 0.25.  Returning  to  the  correlations  of  Figure  30,  we  see 
that  the  zero  crossing  of  the  correlation  for  y * 0.06  6 occurs  at 
tU^/6*  = 24,  which  corresponds  to  a characteristic  frequency  2tt6*/U^t 
= 0.26.  The  correspondence  of  these  frequencies  shows  that  a large- 
scale,  low-frequency  disturbance  coexists  with  small-scale,  high-wave- 
number,  convected  disturbances.  Furthermore,  the  low-frequency 
disturbances  have  a vertical  scale  between  0.22  and  0.44  Inch,  Indicated 
by  the  sign  change  In  the  correlations  of  Figure  31.  Examination  of  the 
flow  pattern  at  position  F (In  Figure  29)  shows  thau  this  distance  from 
the  surface  corresponds  to  the  upper  extremity  of  the  estimated 
separation  region. 

The  further  connection  between  these  measurements  and  the  vortex 
structure  was  made  using  a "visualization"  with  a series  of  cotton  tufts 
which  were  mounted  to  the  strut.  Long  tufts,  2 Inches  In  length, 
moved  In  an  Intermittent  rotary  pattern  which  was  alternated  with  the 
tufts  being  swept  downstream.  Then,  In  an  Instant,  they  resumed  their 
rotary  motion.  Short  tufts,  about  0.5  Inch  In  length,  were  buffeted 
In  an  oscillatory  motion  up  and  dovmstream. 

These  observations  present  a picture  of  an  unsteady  large  vortex 
which  Is  formed  and  then  swept  downstream.  The  core  size  of  this  vortex 
is  probably  on  the  order  Indicated  by  the  closed  curve  in  Figure  29. 

The  process  is  only  weakly  periodic,  occurring  at  the  approximate  fre- 
quency o)6*/U^  = 0.25,  or  wh/U^  = 0.6,  where  h is  the  thickness  of  the 
trailing  edge  as  Indicated  by  velocity  spectra. 


All  the  terms  in  the  equation  have  been  previously  defined,  except  v, 
which  is  the  vertical  (normal)  velocity  fluctuation.  Bldwell  assumed 


that  V vanishes  at  y = 0 and  y » 6.  In  equilibrium  boundary  layers  the 
terms  on  the  right  hand  side  are  negligible,  compared  to  those  on  the 
left  and  so  they  are  neglected.  In  Figure  33  the  mean  momentum  terms 
are  shown  for  the  flow  near  the  edge  at  » 100  feet  per  second.  Up- 
stream of  position  E,  the  dl«fference  between  the  momentum  loss 

(^)  pressure  gradient  terms  agrees  with  the  wall-shear 

coefficient,  determined  by  the  curve-fitting  procedure  outlined  in 
Section  3.  At  position  E,  the  difference  between  the  mean  momentum 
terms  becomes  more  significant  and  indicative  of  the  momentum  lost  from 
the  mean  motion  and  converted  into  the  momentum  of  unsteady  Reynolds 
stress  terms. 


Bldwell,  J.M.,  "Application  of  the  von  Karman  Momentum  Theorem  to 
Turbulent  Boundary  Layers,"  National  Advisory  Committee  for  Aeronautics 
TM  2571  (1951). 


5.2  TURBULENT  FLOW  NEAR  25-DEGREE— 

ROUNDED-BEVELED  EDGE 

The  surveys  of  the  boundary  layer  and  the  static  pressure  distri- 
bution on  the  25-deg — rounded-beveled  edge  are  shown  in  Figure  34. 
Downstream  of  position  F,  the  static  pressure  was  nearly  constant  for 
speeds  U^^^  = 60  and  100  ft/sec;  however,  the  gradient  of  pressure  was 
adverse  downstream  37  Inches  from  the  leading  edge.  The  mean  velocity 
defects  U^-U(y)  (for  U^  * 100  ft/sec)  Increase  markedly  downstream  from 
this  point,  while  the  fluctuating  velocities  Increase  to  an  apparent 
maximum  in  the  vicinity  of  positions  H and  I.  The  balance  of  mean 
momentum  (Figure  35)  shows  that  lost  mean  momentum  is  not  compenstated  by 
the  pressure  gradient  at  locations  downstream  from  position  E-1.  This 
suggests  that  the  boundary  layer  separated  downstream  from  this  point. 

An  examination  of  the  mean  velocity  profiles,  however,  shows  that  the 
normal  gradient  ^ vanishes  somewhere  between  positions  F and  G, 
suggesting  that  separation  occurs  in  that  region.  Apparently,  the 
separation  occurred  near  position  F.  Flow  visualization  by  oil  streaks 
demonstrated  that  flow  separation  occurred  Just  upstream  from  position  F. 
Distinct  flow  reversal  was  not  clearly  disclosed  because  the  very  low 
velocities  near  the  wall  were  insufficient  to  cause  significant  movement 
of  the  oil.  In  drawing  streamlines  in  Figure  34,  position  F has  been 
selected  as  the  separation  point.  Inequality  of  the  mean  momentum 

terms  at  this  position  as  shown  in  Figure  35  is  probably  due  to  the 
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generation  of  Reynolds  stresses  as  Bldwell  has  pointed  out. 

Table  2 gives  the  mean  properties  of  the  boundary  layer  on  the 
rounded  trailing  edge.  Dimensionless  velocity  profiles  are  shown  in 
Appendix  B.  The  favorable  pressure  gradient  upstream  from  position  C 
causes  the  boundary  layer  to  contract  so  that  its  momentum  and  displace- 
ment thicknesses  decrease,  and  the  shape  factor  decreases  to  1.14.  Although 
the  mean  momentum  balance  (Figure  35)  suggests  that  this  is  accompanied  by 
an  Increase  in  C^,  compared  to  the  values  of  local  skin  friction  at 
positions  A and  B.  The  survey  of  the  boundary  layer  was  not  made  with 
close  enough  streamwlse  separations  for  the  calculated  momentum  inte- 


grals to  give  local  values  of  precisely. 


Related  investigations  are  those  of  Schloemer  (1967)  and  of 
28 

Burton  (1973)  in  which  boundary- layer  turbulence  and  wall-pressure 
fluctuations  were  measured  on  surfaces  with  both  favorable  and  adverse 
pressure  gradients.  Although  the  surfaces  on  which  their  measurements 
were  made  were  not  trailing  edges,  the  pressure  gradients  of  their 
measurements  can  be  compared  to  those  of  this  study.  Figure  35  shows 
that  the  favorable  (negative)  gradients  of  the  momentum  balance  in  the 
current  study  are  significantly  greater  in  magnitude  than  those  examined 
earlier,*  while  the  range  of  adverse  (positive)  gradients  examined  in 
the  Burton  measurement  is  representative  of  the  attached  boundary  layer 
on  the  25-deg — rounded-beveled  trailing  edge.  The  favorable  gradient 
of  the  current  dl.  .usslon  is  more  severe  than  the  gradients  of  the 
earlier  two. 

The  frequency  spectra  of  the  longitudinal  velocity  fluctuations  in 
the  boundary  layer  show  the  same  property  of  local  convection  that  has 
been  shown  for  boundary  layers  with  vanishing  pressure  gradients. 

Spectra,  normalized  as  in  Equation  (3),  are  shown  in  Figure  36;  the 
measurements  were  all  obtained  in  the  region  of  adverse  gradient.  At 
position  DE,  the  dimensionless  spectra  at  high  frequencies  are  similar 
to  those  which  were  obtained  on  the  region  of  the  strut  for  which 
dP 

0.  This  similarity  implies  that  the  small  length  scale  is  6*,  which 
is  typical  of  the  high-frequency  turbulence.  Also,  the  dimensionless 
spectra  show  that  low-frequency  motion  which  probably  is  caused  by  large- 
scale  motion  is  reduced  near  the  wall,  i.e.,  near  y ==  0.01A6,  compared 
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Schloemer,  H.H.,  "Effects  of  Pressure  Gradients  on  Turbulent  Wall 
Pressure  Fluctuations,"  Journal  of  the  Acoustical  Society  of  America, 

Vol.  42,  pp.  93-113  (1967). 

28 

Burton,  T.E.,  "Wall  Pressure  Fluctuations  at  Smooth  and  Rough  Walls 
Beneath  a Turbulent  Boundary  Layer  with  Favorable  and  Adverse  Free-Stream 
Pressure  Gradients,"  Massachusetts  Institute  of  Technology  Acoustics  and 
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It 

Since  pressure  coefficients  were  unavailable  in  that  reporting,  the 
approximation  - 0 was  made  in  order  to  calculate  the  momentum  terms. 

47 


11 


to  levels  at  distances  from  the  wall  on  the  order  of  6*.  Near  the  point 

of  separation  at  position  F,  the  dimensionless  spectra  at  different 

heights  above  the  wall  are  separated  Into  two  groups.  For  y/6  > 0.13 

to  0.15  the  turbulence  Is  controlled  by  lower  frequency  (larger  scale) 

disturbances  than  at  points  closer  to  the  wall.  An  alternative  non- 

29 

dlmenslonallzatlon,  suggested  by  Bradshaw  (1967),  replaces  the  dis- 
placement thickness  by  the  distance  from  the  wall.  This  operation  was 
performed  on  the  spectra  of  Figure  36b,  and  It  resulted  In  a collapse  of 
high-frequency  levels  into  two  groups;  the  first  for  y/6  ■ 0.025,  0.033, 
0.296,  and  0.466  and  the  second  for  y/6  ■ 0.127  and  0.152;  see  Figure  37. 

At  the  latter  two  distances  from  the  wall,  the  curvature  of  the  profile 

2 2 

of  mean  velocity  nearly  vanishes,  l.e.,  3 U/3y  = 0;  see  Appendix  B. 

Spectra  which  may  be  nondlmenslonallzed  In  this  manner  have  been 

suggested  by  Bradshaw  to  be  controlled  by  eddies  In  the  near-wall  region 

whose  sizes  are  proportional  to  distance  from  the  wall.  At  the  location 

above  the  wall  where  the  mean  velocity  gradient  Is  maximum,  the  disturbances 

appear  to  be  dominated  by  higher-frequency,  wave  number,  components 

compared  to  those  at  other  locations.  Also  In  contrast  to  the  spectra 

-5/3 

measured  at  position  B,  the  k dependence  occurs  over  a very  restricted 
range  of  wave  numbers.  In  the  region  of  the  separated  boundary  layer, 
position  G,  Figure  36c  shows  that  the  nondimenslonallzatlon  on  6*  gives 
spectra  which  are  generally  similar  throughout  the  boundary  layer.  Also, 
the  spectra  are  similar  to  those  both  at  position  B and  at  position  F 
on  the  knuckle  edge  for  0.8  < U)6*/U(y)  < 6.  For  wave  numbers  outside 
this  range,  the  dimensionless  spectra  at  position  F on  the  knuckle  edge 
are  dissimilar;  however,  they  are  generally  greater  than  those  at 
position  G on  the  round  edge.  Thus,  the  large-scale  motion  on  the 
knuckle  edge  Is  more  Intense  than  that  occurring  In  the  separated  region 
of  the  rounded  edge. 

A review  of  Figures  7,  32,  and  36  shows  that  over  the  range 
1 < (jo6*/U(y)  < 6,  the  dimensionless  spectra  are  generally  similar  and 
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behave  approximately  as  [(i)6*/U(y)]  . This  dependence  on  wave  number 

suggests  the  existence  of  an  Inertial  subrange,  analogous  to  that  which 
exists  in  isotropic  turbulence.  At  higher  frequencies  the  wave  number 
dependence  Increases,  while  at  lower  wave  numbers  the  spectra  occasionally 
show  various  peaks.  In  both  of  these  extreme  ranges,  spectra  are  notably 
dissimilar  at  different  locations  along  the  streamline  and  at  varying 
distances  from  the  wall.  If  we  Interpret  the  dimensionless  wave  number 
as  2tt6*/X^,  where  is  the  streamwise  length  scale  of  convected  eddies, 
the  spectra  of  Figures  32,  33,  and  36  show  that  the  intensity  of  the 
large  turbulent  eddies  for  which  X^  > 2tt6*  are  position-dependent.  The 
lack  of  similarity  of  spectra  at  low  wave  numbers  at  a given  position  on 
either  of  the  trailing  edges  indicates  that  these  larger  eddies  are  not 
convected  at  velocities  proportional  to  U(y).  Within  the  fully  separated 
flow  on  both  edges,  the  dimensionless  spectra  are  similar;  this  would 
suggest  that  6*  is  a spatial  scale  for  the  turbulent  disturbances  at 
some  frequencies. 

6.  PRESSURE  FLUCTUATIONS  ON  NONSINGING  EDGES 
For  both  25-deg — beveled-trailing  edges,  frequency  distribution  and 
spatial  coherence  of  pressures  have  been  measured.  Furthermore, 
coherence  of  pressures  with  turbulence  in  the  locally  separated  flow 
will  be  discussed  to  establish  a mechanism  for  generation  of  the 
pressures. 

6.1  PRESSURES  ON  25-DEGREE— ROUNDED- 
BEVELED  EDGE 

Autospectral  Densities 

The  boundary-layer  flow  on  this  edge  has  just  been  shown  to  be 

both  attached  with  an  adverse  (positive)  pressure  gradient  and  separated. 

The  statistics  of  the  fluctuating  pressures  on  this  edge  are  therefore 

typical  of  those  occurring  in  attached  boundary  layers  with  adverse 

28  27 

pressure  gradients  (Burton  and  Schloemer  ) as  well  as  those  beneath 
fully  separated  flows;  see  Section  6.2. 
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The  region  of  flow  separation  on  this  edge  was  too  small  In  Its 
spatial  extent  to  permit  correlation  measurements.  Thus  a discussion 
of  those  pressure  statistics  will  be  largely  deferred  until  Section  6.2 
so  that  we  will  first  concentrate  on  the  statistics  beneath  the  attached 
portion  of  the  boundary  layer.  In  Figure  38  the  spectral  densities  of 
pressure  fluctuations  have  been  normalized  on  the  free*-stream  dynamic 
head  and  the  local  displacement  thickness.  The  dimensionless  pressure 
spectra  at  positions  A and  B are  identical,  and  they  have  been  shown  to 
be  typical  of  pressures  beneath  boundary  layers  on  smooth  surfaces  in 
zero  pressure  gradient.  Position  C is  in  the  region  of  favorable  gradient, 
and  the  spectrum  levels  are  slightly  lower  than  those  at  position  A. 

For  the  other  measurement  locations  as  the  static  pressure  gradient  in- 
creased, the  pressure  fluctuations  became  more  dominated  by  low-frequency 
disturbances.  At  separation  the  pressure  levels  were  the  highest  at 

oj6*/Uj^  ~ spread  in  dimensionless  spectrum  levels  exhibited  in 

28 

Figure  38  is  a factor  of  5.  Burton  has  shown  spread  of  a factor  of 
10  for  the  range  of  static  pressure  gradients  shown  in  Figure  35.  A 
comparison  of  the  dimensionless  spectra  Just  upstream  of  separation  of 
position  E-1  of  the  current  measurements  and  those  of  earlier  studies 
is  shown  in  Figure  39.  The  spectra  are  designated  according  to  the 
dimensionless  pressure  gradient  parameter  of  Figure  35.  The  spectra  at 
l+H/2 

0 "X-C — ~ notably  dissimilar,  and  the  difference  in  shape  is 

possibly  due  to  the  Influence  of  the  upstream  boundary-layer  history. 

In  comparing  the  current  experiment  with  that  by  Burton,  we  find 

that  in  addition  to  differences  in  N„  at  the  measurement  points,  the 

K 

®0 

momentum  thickness  upstream  from  the  region  of  pressure  gradient 

(position  C)  is  0 = 0.06  inch,  and  C,  = 0.003  upstream  from  position  C, 

28  ^ 

while  for  Burton,  0 = 0.28  inch,  and  * 0.001  at  the  entrance  of  the 
region  of  his  adverse  pressure  gradient.  The  maximum  values  of 
dimensionless  spectra  are  similar  for  both  boundary  layers;  however, 
near  o)6*/U  * 3 the  spectrum  level  on  the  current  trailing  edge  i<5 
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roughly  a factor  of  10  higher  Chan  in  Che  Burton  case.  The  spectra  In 
Figure  38  show  chat  as  Che  fluid  in  the  boundary  layer  of  Che  strut 
moves  Into  Che  region  of  adverse  gradient  (position  DE),  the  pressure 
fluctuations  first  increase  above  Che  upstream,  baseline,  values  at 
zero  pressure  gradient  throughout  the  frequency  range.  Further  down- 
stream the  pressures  at  low  frequencies  continue  Co  Increase  while 
some  reduction  occurs  at  high-dimensionless  frequencies.  Finally  when 
Che  flow  separates  at  position  G,  the  spectrum  is  broad  with  a peak  at 
0)6*/U^j^  * 0.4.  In  the  Burton  case,  a gradually  increasing  adverse 
pressure  gradient  was  also  Imposed  on  Che  boundary  layer  without 
separation;  however,  the  value  of  for  his  incident  flow  was  con- 
siderably less  chan  in  the  current  case.  Thus,  since  pressure  spectrum 

16  18  28 

levels  in  zero  gradient  Increase  with  (Blake,  * Burton,  and 
Section  3 of  this  report)  the  upstream,  baseline,  pressure  levels  were 
notably  lower  in  the  Burton  case  Chan  on  Che  trailing  edge  under  dis- 
cussion. Thus  we  conclude  chat  upstream  history,  which  affects  the 
growth  of  the  incident  boundary  layer  as  well  as  Che  local  Reynolds 
number,  which  is  based  on  local  momentum  thickness,  affects  the  level 
of  local  pressure  fluctuations  on  Che  rounded-trailing  edge  in  Che 
severely  positive  pressure  gradient. 

Figure  39  indicates  Chat  the  maximum  spectrum  level  attained  by 

2 -4 

the  pressure  fluctuations  reaches  the  limit  'I’(w)U^/q^^  6*  = 10  while 
Che  flow  is  attached.  Higher  levels  are  attained  after  separation. 

Spatial  Coherence  and  Convection 

Cross  spectral  densities  of  Che  pressures  at  the  trailing  edge  were 
obtained  with  longitudinally  separated  microphones.  These  spectra  are 
normalized  on  autospectrum  levels  and  are  expressed  as  functions  of 
phase  in  Figure  40.  The  representation  follows  that  of  Section  3 
which  was  used  for  the  boundary  layer  with  zero  pressure  gradient. 

The  normalized  cross-spectrum  represents  a coherence  which  is  not  a 
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unique  function  of  the  phase  a(r  ,co) . A nearly  Identical  behavior  was 

28  ^ 

determined  by  Burton  for  the  smooth  wall  with  adverse  pressure  gradient. 
19 

Corcos,  Interpretation  of  the  longitudinal  phase  of  Equation  (9) 
was  used  again  to  determine  the  convection  velocities  of  Figure  41. 

Data  In  the  upper  figure  were  obtained  about  zero  pressure  gradient  at 
position  B.  They  are  typical  of  measurements  taken  elsewhere.  For 
the  region  of  adverse  gradient,  levels  are  shown  In  the  lower  figure. 

All  separations  are  referred  to  the  lettered  position,  and  they  are  In 
the  direction  of  the  free-stream  flow.  The  convection  velocities  are 
normalized  on  U^;  frequencies  are  expressed  as  In  Figure  38,  except  that 
an  average  of  6*  at  the  measurement  points  Is  used.  A comparison  of 
Figures  38  and  41  shows  that  the  measurements  are  dominated  by  the  flat 
portion  of  the  pressure  spectrum.  The  convection  velocities  for  pressures 

(i)6* 

In  the  frequency  range  rt — >0.1  fall  Into  two  classes,  depending  on 

^00 

whether  the  reference  point  Is  DE  or  E-1.  For  the  teference  position  DE, 

U /U  approaches  0.8;  for  r * 1 Inch,  Figure  40  shows  that  the 
c s 

coherence  exceeds  that  measured  with  E-1,  the  reference.  When  E-1  Is 
the  reference,  the  convection  velocities  approach  0.55  at  the  higher 
dimensionless  frequencies.  We  will  assume  that  the  pressures  are 
generated  by  velocity  disturbances  which  are  located  In  the  boundary  layer 
situated  a distance  y^  from  the  wall  so  that  U(y^.)  “ U^.  An  examination 
of  Figure  37  as  well  as  of  the  dimensionless  data  of  Appendix  B In  light 
of  this  assumption  leads  to  a specification  that  the  pressures  are 
generated  In  the  region  for  which  streamline  velocity  fluctuations  are 
largest.  This  Is  also  the  region  for  which  dU/dy  Is  maximum. 
Interpretation  of  the  cross  spectrum  as  the  coherence  of  a pressure 
disturbance  of  wave  number  leads  to  the  further  speculation  that 

the  disturbance  loses  Its  Identity  as  It  passes  down  the  edge.  The 
Increased  production  of  Reynolds  stresses  as  the  pressure  gradient 
Increases  Is  the  cause  of  this  disorder.  At  frequencies  of  such  a 
nature  that  0)6*/U  < 0.1,  the  convection  velocities  show  two  other 
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distinctly  different  behavior  inodes.  For  the  separation  of  0.5  inch  with 

E-1  the  reference  location,  is  increased  dramatically  with  decreasing 

frequency.  Alternatively,  when  DE  is  the  reference,  the  convection 

velocities  decrease  In  agreement  with  the  Burton  measurements. 

Normalized  cross  spectral  densities  with  laterally  separated 

microphones  are  shown  in  Figure  42  for  separation  distances  ranging  from 

0.31  to  0.62  inch.  The  reference  location  is  position  E-1.  Separations 

are  normalized  on  the  wave  number  (u/U  , where  U is  taken  from  data  in 

c c 

Figure  A1  for  position  E-1  with  r^  ■ 0.5  inch.  The  lateral  cross 

spectral  density  appears  to  be  a universal  function  of  wr^/U^;  it  is 

also  insensitive  to  the  nature  of  the  boundary  layer.  This  is 

Illustrated  by  the  close  similarity  between  the  current  data  and  those 
27 

of  Schloemer  as  well  as  those  for  the  vanishing  pressure  gradient  in 
Figure  42  for  large  separations. 

Overall  spatial  coherence  of  the  pressures  is  shown  by  the  broad- 
band, space-time  correlations  for  streamline  separations  relative  to 
position  E-1.  Figure  43  shows  the  normalized  correlation  function, 
defined  as 


R (;.T)  ■ 

pp  2 2 — 1/2 

[p  (x,t)  p (x+r,t+T)] 


where  in  this  case  x is  the  coordinate  of  position  E-1.  For  both 

positive — the  reference  sensor  upstream  from  the  movable  sensor — and 

negative  separations,  time  delays  for  maximum  correlation  are 

the  same.  Values  of  the  moving-eddy  correlations*  R (r  ,0,T  ) are 

larger  for  positive  separations  than  for  negative  ones.  Convection 

velocities  r /t  U are  approximately  0.7.  The  broadband  correlations 
8 m “ 

are  dominated  by  the  low-frequency  pcessures  which,  according  to 
Figure  41,  are  convected  past  position  E-1  at  relatively  high  velocities. 


* 29 

This  term  is  originally  due  to  Chandiramanl. 
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The  rapid  loss  in  coherence  in  the  moving-eddy  correlation,  which  is  due 
to  the  rapid  boundary  layer  growth,  is  apparently  similar  for  both 
speeds.  Furthermore,  the  moving-eddy  correlation  Is  less  than  that 
previously  measured  in  zero  pressure  gradient. 

To  summarize  the  characteristics  of  spatial  coherence,  we  note  the 
auto  and  cross  spectra  of  the  pressure  fluctuations  have  shown  that 
while  the  Increased  pressure  gradient  at  the  trailing  edge  causes  an 
accelerated  Increase  in  the  boundary  layer  thickness,  the  correlation 
scales  of  the  pressures  relative  to  5*  diminish,  compared  to  situations 
In  zero  pressure  gradient. 

The  lateral  coherence  remains  somewhat  unchanged  by  the  growth  of 

the  boundary  layer.  As  the  boundary  layer  nears  separation  the 

autospectrum  level  Increases  at  low  frequencies,  compared  to  what  it  was 

when  the  pressure  gradient  was  nearly  zero.  Just  upstream  of  boundary- 

layer  separation  at  position  E-1,  the  pressure  field  is  controlled  by 

disturbances  for  which  0)6*/U^  < 1.  Considering  these  disturbances  as 

19  30 

wavelike,  along  the  lines  implied  by  Corcos  and  (1964)  and  examined 

31 

by  Landahl  (1967),  we  Interpret  this  upper  limit  in  frequency  as  a limit 
in  wavelength  for  which  a)6*/U^  ■ ^ This  effectively  places  a 

limit  on  wavelengths  27r/k^  ■ to  greater  than  3.5  6*  with  the  disturbances 
located  at  the  position  in  the  boundary  layer  for  which  dU/dy  is  maximum, 
and  y - 0.7  6*  (0.12  6).  This  is  also  the  location  which  has  been  shown 
to  delineate  between  velocity  disturbances  with  different  dimensionless 
spectral  forms  as  shown  in  Figure  35b. 


Correlation  between  Pressure  and 
Velocity  Disturbances 

In  order  to  establish  a relationship  between  the  fluctuating  wall 
pressures  and  the  local  boundary  layer  turbulence  as  well  as  to  form  a 
basis  for  comparison  between  singing  and  nonsinging  trailing  edges,  a 


Corcos,  G.M.,  "The  Structure  of  the  Turbulent  Pressure  Field  in 
Boundary  Layer  Flows,"  Journal  of  Fluid  Mechanics,  Vol.  18,  pp.  353-378 
(1964). 
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Landahl,  M.T.,  "A  Wave-Guide  Model  for  Turbulent  Shear  Flow,"  Journal 
of  Fluid  Mechanics,  Vol.  29,  pp.  441-459  (1967). 


f 


series  of  pressure-velocity  correlations  was  determined.  Correlations  like 

these  have  been  measured  In  boundary  layers  without  pressure  gradient  by 

32  33 

Wlllmarth  and  Wooldridge  (1963)  and  by  Burton  (1971).  In  the  current 
experiment,  the  space-time  correlations  between  wall  pressure  and  velocity 
disturbances  along  the  streamline  were  determined  for  pressures  at 
position  B,  zero  pressure  gradient;  position  F,  adverse  pressure  gradients 
near  separation;  and  position  G,  fully  separated  boundary  layer.  The 
correlations  are  defined  more  generally  than  Equation  (11)  as 


S 


R (x;;.T)  - P(«.0 

pu  2 2 

[p  (x,t)u  (x+r,t+T)] 
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At  position  B the  correlations  for  r » (0»*'y)  are  practically  the 
same  as  those  measured  In  previous  Investigations  at  positive  values  of 
time  delay;  however,  they  differ  slightly  at  negative  time  delays.  The 
current  measurements  were  made  at  a lower  Reynolds  number  than  In  the 
earlier  work,  and  there  seems  to  be  a consistent  variation,  with  N 

R 

®0 

shown  for  correlations  at  r /6*  • 0.58.  Close  to  the  wall  r * 0.36  6*, 

y y 

the  pressure  is  more  positively  correlated,  compared  to  r^  =•  1.06  6*. 

With  the  larger  separation,  the  pressures  are  more  negatively  correlated 

at  negative  time  delays.  As  Wlllmarth  and  Wooldridge  have  shown, 

correlations  of  this  shape  may  be  Interpreted  as  correlations  of  pressure 

and  streamwlse  velocity  gradient  3u/3x.  Thus,  assuming  that  the  pressure- 

producing  disturbances  are  a function  of  x - U t 

c 


32 

Wlllmarth,  W.W.  and  E.W.  Wooldridge,  "Measurement  of  the  Correlation 
Between  the  Fluctuating  Velocities  and  the  Fluctuating  Pressures  In  a 
Thick  Turbulent  Boundary  Layer,"  Advisory  Group  for  Aerospace  Research  and 
Development  Report  45b  (1963). 

33 

Burton,  T.E.,  "On  the  Gei.eratlon  of  Wall  Pressure  Fluctuations  for 
Turbulent  Boundary  Layers  over  Rough  Walls,"  Massachusetts  Institute  of 
Technology  Acoustics  and  V'ibiatlons  Laboratory  Report  70208-4  (1971). 


55 


^ OIV 

3—  R (x;r,x)  - -1/U  -5—^ 
3t  pu  c dr 

X 


-1/U  R (x;r,T) 
c pu 

X 


It  can  be  seen  by  Inspection  of  Figure  44  that  R^^  will  have  a maximum 
near  T = 0,  u ■ 3u/3x.  ^ 

X 
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According  to  the  Kraichnan  (1956)  analysis,  the  fluctuating  wall 
pressures  beneath  homogeneous  boundary  layers  are  gererated  by  a source 

field  of  volume  density  p ^ , which  is  continuously  distributed 

throughout  the  boundary  layer.  This  "turbulence  mean  shear"  interaction 

32  30 

has  been  argued  by  Wlllmarth  and  Wooldridge  and  Corcos  to  be  responsible 
for  a strong  correlation  at  T * 0 between  pressure  and  the  streamline 

gradient  of  the  normal  velocity  . This  correlation  they  showed  to  be 

strongly  positive,  and  it  accompanies  a negative  correlation  maximum  with 

3u/3x.  The  space-time  correlations  of  Figure  44  ere  consistent  with  the 

32 

Wlllmarth  and  Wooldridge  results  for  < 6*.  For  r^  ■ 1.06  6*,  how- 
ever, the  shape  of  the  correlation  function  suggests  that  an  alternative 
pressure-producing  mechanism  exists  for  eddies  in  the  outer,  mixing, 
region  of  the  boundary  layer.  It  will  be  shown  that  this  function  shape 
will  be  typical  of  the  separated  boundary  layer  and  so  will  be  further 
discussed  at  the  end  of  this  section. 


Pressure-velocity  correlations,  which  were  obtained  at  position  F 
for  r ■ (0,r  ) are  shown  in  Figure  45  for  * 60  and  100  feet  per  second. 
At  position  F the  maximum  intensity  of  the  turbulence  occurs  at  y ■ 0.54  6* 
for  “ 100  ft/sec,  and  y * 0.52  6*  for  * 60  feet  per  second. 

The  general  functional  dependence  on  T as  well  as  the  maximum  correlations 
is  similar  to  position  B.  At  a corresponding  location  in  the  boundary 
layer  r^  = 0.5  6*,  the  maximum  coherence  is  reduced  at  position  F, 
compared  to  that  at  position  B.  The  trend  of  more  negative  correlations 
with  increasing  r^  is  also  prevalent  at  this  position. 


Kraichnan,  R.H.,  "Pressure  Fluctuations  in  Turbulent  Flow  over  a Flat 
Plate,"  Journal  of  the  Acoustical  Society  of  America,  Vol.  28,  p.  378  (1956), 


Correlations  obtained  above  and  downstream  from  position  G at 

U = 60  ft/sec  are  shown  in  Figure  46.  In  these  cases  r = (r  ,r 
00  ° X y 

=0.1  in.),  where  r^  is  variable.  Lettered  locations  a through  d 
correspond  to  values  of  r^,  which  are  Indicated  in  Figure  37.  The  locus 
of  correlations  Is  just  outside  the  region  of  separation.  These  corre- 
lation functions  are  entirely  negative;  the  level  of  correlation  is  high 
for  positions  a through  d,  and  it  decreases  at  position  e.  The  temporal 
width  of  the  correlation  peak  At,  measured  for  coherence  to  decrease  to 

OO 

one-half  Its  maximum  value,  Is  ~ 8,  where  6*  is  evaluated  at 

ti)6* 

position  G.  This  corresponds  to  a dimensionless  frequency  — = 0.8, 

^00 

which  is  comparable  to  the  bandwidth  of  the  autospectrum  of  fluctuating 

pressure  at  position  G;  see  Figure  38.  On  the  basis  of  this  similarity 

in  frequency,  it  appears  that  the  pressure  generated  beneath  the 

separated  boundary  layer  is  caused  by  a large-scale  velocity  field,  having 

a coherent  influence  that  persists  slightly  downstream  from  the  trailing 

edge.  Also  it  appears  from  the  functional  form  of  the  correlation  that 

the  mechanism  of  pressure  generation  by  the  turbulence  is  different  for 

the  separated  boundary  layer  than  it  is  for  the  attached  boundary  layer. 

The  magnitude  of  the  absolute  value  of  the  correlation  is  also  greater 

for  the  separated  than  it  is  for  the  attached  boundary  layer. 

The  outer  flow  at  position  B is  typical  of  that  existing  in  fully 

developed  boundary  layers,  and  it  is  characterized  by  a nearly  vanishing 

35 

mean  shear;  it  has  been  called  wakelike  by  Coles  (1956).  Although 
this  wake  analogy  comes  about  because  the  mean  velocity  defect  has  certain 
functional  similarities  with  those  of  fully  developed  wakes,  a wake 
analogy  may  also  apply  to  the  mechanism  of  turbulent  pressure  generation. 

A large-scale  turbulent  eddy  structure  is  apparently  convected  at  a 
velocity  nearly  equal  to  over  the  surface  and  trailing  edge.  The 
streamwlse  components  of  the  velocity  field  due  to  the  eddies  are 
negatively  correlated  to  the  pressure  on  surface.  In  Section  4 it  was 


^Coles,  D.,  "The  Law  of  the  Wake  in  the  Turbulent  Boundary  Layer," 
Journal  of  Fluid  Mechanics,  Vol.  1,  pp.  191-226  (1956). 


shown  that  the  surface  pressures  induced  by  the  wake  vorticity  were 
negatively  correlated  with  the  streamwise  component  of  velocity;  see 
Figures  16,  22,  and  23.  For  the  edges  considered  thus  far,  pressures  at 
singing  and  nonsinging  edges  are  distinguished  by  the  degree  of  coherence 
in  the  separation  of  the  boundary  layer.  This  is  reflected  in  the  band- 
width of  the  pressure  autospectrum  as^  well  as  in  the  coherence,  and  the 
downstream  perslstance  of  the  coherence,  of  the  pressure  wi.'-h  the  velocity 
fluctuations.  This  coherence  is,  in  turn,  dependent  on  the  coherence  of 
the  velocity  fluctuations  in  the  separated  boundary  layer  and  wake. 

The  shapes  of  the  pressure-velocity  correlations  at  position  B for 
T < 0 may  be  due  to  the  Influence  of  a large-scale  eddy  structure, 
generated  by  the  laminar  separation  at  the  leading  edge  and  convected 
downstream.  The  mechanism  of  the  pressure  production  is  the  same  as 
that  implied  by  the  correlations  described  prevlouslyt 

6.2  PRESSURES  ON  25-DEGREE— KNUCKLE- 
BEVELED  EDGE 

The  flow  field  on  this  trailing  edge  has  been  shown  in  Section  5 to 
be  partially  representative  of  the  separated  boundary  layer  downstream 
from  position  G on  the  25-deg — rounded-beveled  edge.  The  scale  and 
intensity  of  the  boundary  layer  and  its  velocity  disturbances  was  shown 
to  be  larger  on  the  knuckle  edge.  This  section  is  therefore  an 
examination  of  the  statistical  structure  of  the  pressure  field  beneath 
a separated  boundary  layer.  The  scope  of  the  measurements  is  the  same 
as  that  for  the  rounded  edge. 

Autospectral  Densities 

The  dimensional  autospectra  G(f)  are  shown  in  Figure  47;  they  were 
obtained  at  the  lettered  positions  referenced  in  Figure  29.  At 
U^  = 50  ft/sec,  the  autospectra  at  positions  A through  D are  similar 
within  a factor  of  two;  this  is  expected  because  the  boundary-layer — 
displacement  thicknesses  are  comparable  for  all  these  positions;  see 
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Tables  1 and  2.  At  position  E,  on  the  downstif^eam  side  of  the  knuckle, 

the  pressures  are  more  dominated  by  low-frequency  disturbances.  At 

positions  E-1  and  F,  the  boundary  layer  is  fully  separated,  and  the 

pressures  at  low  frequencies  dominate;  they  are  comparable  for  both 

locations.  At  ■ 100  ft/sec,  the  pressures  are  again  similar  at 

positions  A through  C and  attain  Intermediate  levels  at  positions  D and 

E.  At  positions  E-1,  F,  and  G the  autospectra  are  nearly  Identical  at 

low  frequencies  with  some  spread  for  f > 1 kilohertz.  Beneath  the 

separated  boundary  layer,  the  low-frequency,  pressure  spectrum  levels 

3 

are  factors  of  30  and  10  higher  than  the  pressures  beneath  the 
attached  flow  at  speeds  of  50  and  100  feet  per  second,  respectively. 

For  comparison  of  these  pressures  to  those  of  the  other  flow  situ- 
ations, the  autospectra  are  nondimensloiialized  on  the  outer  flow 

variables  and  6*  in  Figure  A8.  It  can  be  seen  that  the  dimensionless 

2 

broadband,  low-frequency  pressures  are  a factor  of  nearly  10  higher 
than  those  at  position  G on  the  25-deg — rounded-beveled  edge.  The 

(i)6* 

spectrum  levels  in  the  separation  region  diminish  with  rj — > 1 in 

^OO 

similarity  to  those  on  the  25-deg  rounded  edge  at  position  G.  A notable 

feature  of  the  dimensionless  autospectra  beneath  the  separated  boundary 

layer  on  the  knuckle  edge  is  that  the  levels  do  not  scale  on  the  free- 

stream  dynamic  head,  while  the  spectrum  shapes  are  quite  similar  when 

expressed  in  terms  of  u)6*/U^.  That  the  spectrum  levels  are  not  dependent 
2 

on  q“  is  surprising  in  light  of  the  fact  that  the  normalized  velocity 
~2  1/2 

disturbances  u /U^  are  similar  for  both  speeds;  see  Appendix  B. 


Correlations  of  Wall  Pressures 

Space-time  correlations  of  wall  pressures  on  the  knuckle  trailing 
edge  are  shown  in  Figure  A9.  The  definition  of  the  correlation  function 
has  been  given  by  Equation  (2A) . As  before,  a separation  vector  pointed 
upstream  is  negative;  for  the  correlations  of  Figure  A9,  position  F is 
the  referenced  location.  The  coherence  between  pressures  at  positions 


t 


I 


E-1  and  F (r^  = -1  in.)  Is  nearly  0.5  at  speeds  = 50  and  100  feet 
per  second.  The  mean  convection  velocity  is 


U 

c 


s m “ 
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which  is  considerably  less  than  that  which  has  been  calculated  for 
pressures  beneath  fully  attached  boundary  layers.  The  streamwise  and 
lateral  spatial  correlations  R (r. ,0),  where  r.  = r or  r , are  also 
shown  in  Figure  49.  Overall,  the  correlations  show  the  existence  of  a 
large-scale  convected  pressure  field  that  seems  to  be  laterally  corre- 
lated over  two  trailing-edge  thicknesses  while  Instantaneously  extending 
only  a 1/3-trailing-edge  thickness  upstream.  This  spatial  extent 
roughly  corresponds  to  the  estimated  upstream  boundary  of  the  separation 
region;  see  Figure  29.  A correlation  between  pressures  at  positions  F 
and  G shows  convection  in  two  directions.  As  shown  in  Figure  50  two 
contributions,  correlated  nearly  equally,  appear  to  be  convected  at 
U = 0.35  U . The  functional  form  of  the  correlation  Indlca. ss  that  the 
correlation 


p(Xp,t)  1^  (Xg,t+T) 

is  a maximum  at  T * 0.  This  means  that  the  pressure  at  position  F is 
positively  correlated  with  the  temporal  rate  of  change  of  pressure  at 
position  G.  Comparison  of  the  correlation  magnitude  between  pressures  at 
F and  G with  that  of  pressures  at  F and  E-1  (r^  = + 1 inch)  shows 
comparatively  less  correlation  between  pressures  at  F and  G.  In 
Section  5-1,  the  unsteady  character  of  the  velocity  fluctuations  near 
the  edge  was  described  as  due  to  an  unsteady  large  vortex,  alternately 
formed  and  swept  downstream  by  formation  and  ejection  of  the  vortex. 
Positive  correlation  between  p(Xp,t)  and  p(x^,t)  suggests  that  pressures 
at  F and  G are  related  to  the  formation  of  the  vortex  and  its  acceleration 
proceeding  ejection.  They  do  not  appear  to  be  caused  by  a simple  vortex 
convection. 
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Pressure  and  Velocity  Correlations 

Further  indications  of  the  pressure-producing  mechanism  of  the 
separated  boundary  layer  are  given  by  correlations  defined  by 
Equation  (25)  with  r = (0,r  ).  The  space  and  time  correlations  are 
entirely  negative  as  shown  in  Figure  51.  The  lettered  locations  denote 
the  anemometer-probe  positions  designated  in  Figure  29.  At  r = 0.596*, 
the  correlation  nearly  vanishes,  while  it  is  negative  at  r = 0.0346*  and 
r > 1.186*.  The  increase  in  the  absolute  value  of  the  correlation  with 

y - 

increasing  r^  is  somewhat  a result  of  normalization  on  local  turbulence 
intensity.  This  latter  factor  also  decreases  with  large  r^;  see 
Figure  29  and  Appendix  B.  ’ The  maximum  value  of  the  pressure-velocity 
correlation  occurs  for  r^  = 1.186*.  That  the  pressure  is  partially 
generated  by  a convectlng  vortex  field  is  seen  by  reference  to  the 
correlation  in  Figure  52.  The  locus  of  points  g through  k roughly 
corresponds  to  local  maxima  in  turbulence  intensity  as  shown  in  Figure  29. 
It  is  seen  that  while  intensity  decreases  with  increasing  r^,  coherence 
increases  slightly.  In  comparing  the  correlations  in  Figures  46  and  51, 
we  see  that  the  pressures  on  the  25-deg  knuckle  edge  are  much  more 
coherently  generated  by  the  near  wake  than  those  on  the  25-deg — rounded- 
beveled  edge.  Also,  the  pressures  are  negatively  correlated  with  the 
streamwlse  velocity  disturbances  in  similarity  with  correlations  on 
singing  edges.  Figures  16,  22,  and  23.  The  spatial  correlations  of 
pressure  velocity,  Rp^(r,0)  show  that  = 0 at  point  i,  which  is  just 
above  the  tip  of  the  edge.  Here,  the  correlations  change  from  positive 
to  negative  values,  which  is  a behavior  similar  to  that  shown  for  the 
45-deg — rounded-beveled  edge  in  Figures  22  and  23.  The  magnitude  of  the 
spatial  coherence  is  comparable,  i.e.,  R = -0.02.  On  the  knuckle 
edge  the  disturbances  (Figure  31)  are  generally  broadband,  showing  a 
small  low-frequency  peak  at  ~ 0.15  to  0.2 — ^ = 0.5,  where  h is  the 

'^OO  ^00 

thickness  of  the  strut — which  seems  to  be  more  pronounced  at  position  H 
for  y = 0 than  at  position  F. 


The  sets  of  correlations  in  Figures  49  through  52  show  the  existence 

of  a pressure  field  which  Is  partially  generated  by  a wake  that  Is  less 

structured,  compared  to  that  causing  the  periodic  pressures  on  singing 

trailing  edges.  The  disturbances  do  not  grow  as  they  move  downstream  In 

the  wake  as  do  those  In  the  flow  region  of  the  singing  edges;  however, 

they  are  generated  by  the  local  vortlclty  caused  by  the  separation  through 

a mechanism  of  pressure  generation  that  appears  to  be  common  to  both 

singing  and  nonsinging  trailing  edges.  Although  the  magnitudes  of 

pressures  which  are  generated  beneath  separated  boundary  layers  on  both 

the  knuckle  and  the  rounded  trailing  edges  are  quite  different,  to  some 

degree  their  mechanisms  of  generation  seem  to  be  similar  as  can  be 

seen  by  reference  to  the  pressure-velocity  correlations  of  Figures  46  and 

51,  53,  and  54.  For  both  edges  the  pressure-velocity  correlations  In 

Figures  53  and  54  are  shown  for  corresponding  locations  on  the  edges 

relative  to  estimated  regions  of  separation.  The  correlations  are  In 

the  form  R (r  ,0).  For  both  edges,  R (r  ,0)  corresponds  In  shape  to 
y AU/U^  y ~2  1/2 

the  profile  of  mean  shear  -r — as  well  as  to  profiles  of  u /U  versus 


the  profile  of  mean  shear  as  well  as  to  profiles  of  u"  /U^  versus 

y/6*.  The  vertical  extent  of  the  region  of  maximum  correlation  appears 

to  be  measured  by  the  reciprocal  of  the  shape  factor  6*/6.  Finally 

(R  ) for  the  knuckle  edge  Is  three  times  larger  than  for  the  25-deg — 

' pu  max  6 0 6 

rounded-beveled  edge.  This  similarity  In  with  the  velocity 

gradient  Indicates  that  the  source  distributions  above  the  two  edges  are 
as  similar  as  are  the  velocity  profiles.  The  time  behavior  of  each  of 
the  local  correlation  functions  suggests  the  simplified  Interpretation 
that  pressure-producing  disturbances  are  convected  at  a constant  velocity 
by  a fixed  point,  and  a locally  positive  (or  Increased)  velocity  dis- 
turbance Is  associated  with  a negative  pressure  disturbance.  This 
Interpretation  Is  drawn  from  the  theory  of  potential  flow.  However  this 
Is  only  a contribution  to  the  total  pressure.  Other  contributions  are 


generated  by  the  disturbances  ejected  into  the  wake  so  that  the  overall 
level  of  the  pressure  fluctuations  is  generated  by  the  coherent  effect 
both  of  Che  convected  wake  field  and  of  the  local  shear  layer. 


Cross-Spectra  of  Pressure  and 
Velocity  Disturbances 


To  clarify  some  of  the  trends  in  correlations,  cross  spectral 
densities  of  pressures  at  two  streamline  locations  and  cross  spectra 
between  pressure  and  velocity  fluctuations  were  made  to  show  contri- 
butions due  to  various  parts  of  the  near-  and  far-wake  flow  layer. 
Magnitudes  and  phases  of  normallzed-pressure  cross  spectral  densities  are 


shown  in  Figure  55  for  a variety  of  streamline  separations.  The  phase 
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has  been  given  Corcos  Interpretation  as  stated  in  Equation  (9a). 


Cross  spectra  (Figures  55a  through  55c)  show  that  two  frequency  ranges 
of  correlation  exist.  Overall,  from  positions  E to  F,  the  coherence 
is  greatest  near  40  Hz  (oih/Ug^  ■ 0.43),  which  roughly  corresponds  to  the 
frequency  of  the  broad  peaks  in  the  velocity  spectra  of  Figure  32. 
Between  positions  E-1  and  F,  the  coherence  is  greatest  near  200  Hz 
(wh/U^  ” 2.1).  This  frequency  is  a limit  below  which  the  autospectrum 
of  pressure  maintains  its  asymptotic  low-frequency  value.  Locally 
(Figure  55d)  the  pressures  at  position  F are  well  correlated  throughout 
the  frequency  range.  The  disturbances  are  convected  between  positions 
E and  F at  speeds  U = 0.32  and  between  E and  E-1  at  U = 0.35  U ; 


however,  between  E-1  and  F,  = 0.48  U^.  Near  position  F,  =0.4  to 


0.6  Ugg,  depending  on  frequency.  These  velocities  are  consistent  with 
those  indicated  by  the  broadband  correlations.  Convection  velocities 
of  approximately  0.3  are  somewhat  less  than  the  mean  velocities  in  the 
boundary  layer  at  locations  of  the  maximum  streamline  velocity  fluctuations 
which  are  from  0.5  to  0.6  U^.  These  low-frequency  pressures  are  therefore 
convected  at  speeds  more  representative  of  the  average  mean  velocity  of 
the  large-scale  eddy  structure  beneath  the  separated  shear  layer.  The 
more  rapidly  convected  high-frequency  disturbances  have  apparent  sources 
at  the  location  of  maximum  streamline  velocity  fluctuations. 


The  cross  spectra  between  pressures  at  F and  the  local  velocities 
at  the  points  g,  h,  and  1 In  Figure  29  also  show  two  distinctly  different 


frequency  ranges  of  coherence.  In  Figure  56,  for  “ 0*5  inch  and 


r^  >•  2 Inches,  point  1,  the  dependence  of  phase  a versus  frequency,  has 


two  slopes;  the  magnitude  of  the  cross  spectrum  has  two  peaks.  The 


high-frequency  peak  is  convected  at  = 0.9  U^.  At  lower  frequencies. 


the  slope  da/3b0  is  negative;  which  implies  the  existence  of  a negative 


convection.  This  interpretation  follows  a special  definition  of 
convection  velocity  as  * r • Strictly,  is  analogous  to  a 


group  velocity  and  becomes  identical  to  (Equation  (9a))  when  is 


Independent  of  frequency.  For  r^  = 1 inch,  point  h,  the  phase  is  small 


and  nearly  Independent  of  frequency;  for  r^  * 0,  point  g,  the  phase  is 


negative  at  all  frequencies.  In  the  last  case  the  resolvent  separation 
vector,  which  is  vertical,  has  a component  in  the  plane  of  the  edge 


which  is  directed  upstream,  r^  ■ -0.21  inch.  The  convection  velocity  is. 


therefore,  small  and  positive  as  designated  in  Figure  56.  This  is 
representative  of  the  low  speed  of  inner  entrained  fluid.  That  the  phase 
is  small  for  point  g suggests  that  the  pressure-producing  disturbances 
in  the  flow  pass  points  F and  G at  the  same  time.  The  similarity  between 
convection  velocities  in  the  cross  spectral  densities  of  Figures  55  and 
56  show  that  low-frequency  pressures  are  generated  by  the  low-speed 
large  eddy  flow,  while  the  higher  frequency  pressures  f > 100  Hz  are 
basically  generated  by  the  smaller  scale  eddies  in  the  shear  layer. 


Comparisons  of  Singing  and  Nonsinging 
Edges 


Finally,  the  degree  of  dependence  of  the  fluctuating  pressures  on 
the  eddy  structure  in  the  wake  has  been  further  indicated  by  the  cross 
spectra  of  pressures  on  opposite  sides  of  the  knuckle  trailing  edge. 
Sensors  were  positioned  at  E-1  (Figure  29)  on  opposite  sides  of  the 
trailing  edge,  and  the  cross  spectra,  measured  in  12.5-Hz  frequency  bands 
are  shown  in  Figure  57.  The  high  level  of  coherence  for  low  frequencies 


further  supports  the  hypothesis  that  these  pressures  are  generated  by  the 
large-scale  disordered  wake  field.  The  frequency  has  been  made 


dimensionless  on  h and  U_.  This  thickness  Is  selected  here  as  a 


convenient  scale  which  Is  representative  of  the  knuckle-edge — shear-layer 
separation,  and  It  Is  used  to  draw  a comparison  with  the  singing  spectra 


In  Figures  17  and  21.  For  frequencies 


> 1 and  t*)h/U  > 1 (taking 


y^  - h,  = U^)  all  the  pressure  spectrum  levels  decrease.  It  Is  to  be 


noted  that  the  narrowband  dimensionless  pressures  In  Figures  15,  21,  and 
57  may  all  be  compared  directly.  It  Is  Interesting  that  even  though  the 
25-deg  knuckle  edge  did  not  generate  a periodic  vortex  street  In  the 
speed  range  of  the  measurements.  It  did  generate  high-level  pressure 
fluctuations  In  a frequency  range  not  far  removed  from  the  frequency  at 
which  singing  would  have  occurred — wh/U^  of  order  unity. 

In  the  case  of  the  25-deg  rounded  edge,  the  pressures  at  position  G 
are  considerably  less  than  those  which  are  generated  by  other  edges; 


nonetheless,  they  are  restricted  to  frequencies  toy^/U^  < 1,  where  y^  Is 


the  distance  between  maxima  In  turbulent  velocity  Intensities  shown  In 


Figure  35,  y^  - 0.5  Inch.  A complete  series  of  measurements  was  not 


made  on  the  25-deg  rounded  edge  at  U ■ 200  ft/sec,  however,  at  this 


speed,  the  pressures  at  position  F began  to  display  a weak  tone  at 
f - 450  hertz.  Assuming  that  y^  = 0.5  Inch — suggested  by  the  data  at 


positions  H and  I In  Figure  35 — and  that  this  frequency  Is 


equivalent  to  ooy^/U^  = 1.2.  Weakly  periodic  velocity  disturbances  w^ere 


also  measured  In  the  wake  at  this  speed. 


7.  SUMMARY:  TRAILING  EDGES  AND  UNSTEADY  PRESSURES 


We  have  seen  that  all  the  edges  considered  In  this  report  generate 
high-level — f luctuatlng-surface  pressures,  regardless  of  whether  or  not 
periodic  disturbances  are  formed  In  their  wakes.  It  has  also  been  shown 
that  the  periodic  disturbances  occur  with  Reynolds-number-dependent 
strengths  on  the  blunt  edge  and  on  the  45-  and  25-deg — rounded-beveled 
edges.  The  flow  near  the  knuckle  edge  was  not  examined  over  a large 
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enough  range  of  speeds  to  determine  if  periodic  pressures  would  occur  at 
high  Reynolds  number.  When  periodic  disturbances  were  not  generated,  the 
fluctuating  pressures  on  the  edge  were  distributed  over  a frequency  range 
bounded  by  the  upper  limit  of  frequency  U)y^/U  = 1.  The  spanwlse 
correlation  length  of  the  pressure  fluid  on  the  knuckle  edge  was  one-fifth 
that  of  the  periodic  pressure  field  on  the  blunt  edge.  Both  the  low- 
frequency  random  pressures  on  the  one  hand  and  the  periodic  pressures 
due  to  coherent  vortex  streets  on  the  other  hand  appear  to  be  generated 
by  large-scale  vortlclty  through  a common  mechanism.  Higher  frequency 
pressures  In  the  random  cases  appear  to  be  due  to  a smaller-scale 
convected  vortex  structure,  which  Is  located  at  the  upper  boundary  of 
the  region  of  separation.  Contributions  to  the  random  velocity  disturbances 
In  the  separated  and  nearly  separated  flow  on  nonsinging  edges  which  are 
locally  convected  appear  to  be  due  to  a vortex  structure  whose  length 
scale  Is  governed  by  the  local  displacement  thickness.  The  other  con- 
tributions to  the  velocities  which  are  not  locally  convected  and  which 
appear  to  be  generated  by  the  near-wake  structure  occur  at  frequencies 
dependent  on  the  shear-layer  separation  In  the  near  wake. 

A comparison  of  the  turbulence  Intensities  and  length  scales  of  all 
edges  shows  that  no  correspondence  exists  between  the  degree  of 
periodicity  of  the  surface  pressures  and  the  value  of  the  maximum 


broadband  Intensity  of  velocity  fluctuations  at  the  edge  extremity. 


“2  1/2 

These  Intensities  cover  the  range  0.16  < u /U  < 0.23.  Periodic 

- o - 


pressures  have  been  shown  by  the  experiments  to  require  the  growth  of 


coherent  wavelike  disturbances  in  the  wake,  and  this  growth  Is  therefore 

-T  1/2 


apparently  unrelated  to  u 


/U^.  This  evidence  has  been  supported  by 


the  stability  calculations  of  Appendix  C which  show  the  amplification 
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factors  of  disturbances  to  be  dependent  on  Reynolds  number  and  d U/dy  . 


The  remainder  of  this  section  will  be  devoted  to  the  practical 
questions  of 

1.  Establishing  whether  or  not  a given  trailing  edge  will  generate 
a periodic  vortex  street 


2.  Estimating  the  frequencies  of  pressures  that  are  generated 

3.  Estimating  fluctuating  pressure  levels  and  length  scales. 

Although  certain  dimensionless  factors  have  been  shown  in  the  preceeding 
sections  to  describe  the  measurements  quite  precisely,  these  quantities 
are  probably  not  very  useful  for  estimations  because  they  Involve 
parameters  which  are  associated  with  flow  details  rather  than  with 
overall  geometric  details.  In  an  estimation  process  the  quantities  are 
not  known  £ priori,  and  they  are,  therefore,  of  little  practical  value. 

The  results  of  the  measurements  will  be  summarized  and  approximated 
below  so  that  some  simple  formulas  may  be  given. 

Minimum  Reynolds  Number  for  Periodic 
Wake  Generation 

Periodic  disturbances  were  generated  on  the  blunt  and  both  of  the 

rounded-beveled  edges  above  certain  onset  speeds.  Using  a definition  of 

Reynolds  number  based  on  the  wake  parameter  y.,  these  speeds  corresponded 
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to  yr/v  ■ 5.2  X 10  for  the  25-deg  edge,  3.3  x 10  for  the  45-deg 
* 3 

edge,  and  4.1  x 10  for  the  blunt  edge.  Length  scales  are  y^  « 0.25  h 
and  0.6  h for  the  25-  and  45-deg  round  edge,  respectively,  and  0.75  h for 
the  blunt  edge.  Although  the  onset  Reynolds  ntmibers  are  not  the  same 
for  all  edges,  there  appears  to  be  a consistent  Increase  in  onset 
Reynolds  number  as  the  apex  angle  of  the  edge  is  decreased. 

Vortex  Shedding  Frequencies  • 

Although  a universal  Strouhal  number  based  on  a geometric  length  of 
the  trailing  edge  has  not  been  Isolated  in  this  study,  an  approximate 
dimensionless  frequency  is  f^  yj/^oo  “ 1/2it.  The  length  y^  for  the  blunt 
and  rounded-beveled  edges  has  been  given  above.  Other  dimensionless  forms 
have  been  proposed  by  other  authors  for  specific  tralllng-edge  forms,  and 
they  may  give  frequencies  more  precisely  in  those  cases;  however,  none 
appear  to  have  practical  and  universal  application.  When  shedding  does 
not  occur,  the  random  pressures  are  bounded  by  the  upper  limit  in  freqr.aucy 
~ Here,  values  of  y^  are  the  same  as  those  given  previously; 

for  the  25-deg  knuckle-beveled  edge,  it  is  approximately  0.6  to  0.8  h. 


- A. 


Fluctuating  Pressure  Levels 

The  surface  pressure  levels  exerted  on  the  struts  apparently  are 
not  governed  by  any  known  simple  scaling  law.  In  general,  at  a given 
speed  the  level  of  the  periodic  pressures  on  the  rounded-beveled  edges 


decreases  with  a reduction  In  apex  angle.  When  random  pressures  were 
generated  on  the  rounded  edges,  the  broadband  mean-square  levels  were 

-t^  2 

seen  to  be  comparable  to  10  q^.  For  the  25-deg — knuckle-beveled  edge, 

the  mean-square  levels  of  the  random  pressures  attained  substantially 

-2  2 

higher  levels,  approaching  10  q^. 
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Figure  2 - Ratio  of  Response  of  Probe  Tubes  to  Response  of  a 

1/2-Inch  Microphone 
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Figure  6 - Mean  and  Turbulent  Velocity  Profiles  for  Positions  4 
Open  Points,  and  B,  Closed  Points,  on  Strut  with  25-Degree — 

Rounded-Trailing  Edge 


Figure  7 - Spectra  of  Longitudinal  Turbulent  Velocity  at  Position 
on  Knuckle  Edge  at  U * 100  Feet  per  Second 
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Figure  9 - Narrowband  Mean-Square  Pressures,  Measured  in  62.5-Hertz  Bands 
at  Forward  Chordwlse  Locations  of  Knuckle-Trailing  Edge 
at  U ■ 50  Feet  per  Second 
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Figure  12  - Normalized  Magnitude  of  Longitudinal  Cross  Spectral  Density 

at  Position  B 
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Figure  16  - Pressure  Velocity  Correlations  behind  Blunt  Edge  at 
« 100  Feet  per  Second;  r^  Is  Measured  Relative  to  Edge 

(In  diagram  (d)  open  points  were  obtained  from  R ; 
cross  points  from  R .) 
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Figure  18  - Lateral  Spatial  Correlation  of  Fluctuating  Traillng-Edge 
Pressure,  Measured  at  100  Feet  per  Second 


Figure  21  - Mean-Square  Pressures,  Measured  on  the  45-Degree  Rounded  Edge 
in  12.5-Hertz  Bands  at  Positions  G and  6,  Flagged  Points; 
for  = 50  Feet  per  Second,  Open  Points,  Awy^/U^  * 0.09; 

and  for  100  Feet  per  Second,  Closed  Points, 


)x  10 
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DISTANCE  FROM  TRAILING  EDGE  TIP,  r 

*/IN. 

Figure  23  - Pressure-Velocity  Correlations  behind  45-Degree 
Rounded-Beveled  Edge  at  = 100  Feet  per  Second; 

Distances  r^  are  Measured  from  Tip  of  Edge 

(Points  marked  0 along  y * y^t  A along  y = y^^) 
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'ortex  Geometry  in  Near  Wake,  Corresponding  to  Pressure' 
Velocity  Correlations  for  Upper  Shear  Layer 
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Figure  28  - Universal  Strouhal  Number  of  Bearman  (1967)  versus  U /U 

s 

for  Blunt  and  45-Degree — Rounded-Trailing  Edges 
(A  45-degree  round,  100  feet  per  second,  0 blunt, 

100  feet  per  second; — representative 

body  shapes.) 


Figure  31  - Space-Time  Velocity  Correlations  above  Position  F on 
25-Degree — Knuckle  Edge,  Using  Vertical  Separations  of  Probes 
(U  “ 50  feet  per  second;  6 = 1,69  inch.) 


Figure  32  - Spectra  of  Longitudinal  Turbulent  Velocity  at 
Positions  F and  H on  Knuckle  Edge 
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Figure  32a  - Position 
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Momentum  Balance  for  25-Degree  Knuckle  Edge 
at  U * 100  Feet  per  Second 
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Figure  36  - Autosp  Ta  of  Turbulent  Velocity  Fluctuating 
over  25-Degree — Rounded-Beveled  Edge 
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Figure  37  - Autospectra  at  Position  F on  25-Degree — Rounded 
Beveled  Edge  Nondimensionalized  on  Distance  from  Wall 
(Data  are  taken  from  Figure  36b.) 
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Comparison  of  Various  Measurements  of  Autospectra  of 
ill  Pressures  In  Adverse  Pressure  Gradients 
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Figure  42  - Normalized  Lateral  Cross  Spectral  Density  Magnitudes 
of  Fluctuating  Surface  Pressure  at  Position  El  on 
25-Degree — Rounded-Beveled  Edge 
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Figure  43  - Broadband  Space- Time  Correlations  of  Pressures  at 
Position  El  on  25-Degree — Rounded-Trailing  Edge 
(6*  = 0.12  inch  selected  at  position  El  for  * 100  feet  per  second) 

( 60  feet  per  second,  100  feet  per  second; 

autocorrelations:  0 60  feet  per  second,  0100  feet  per  second) 
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Figure  44  - Pressure-Velocity  Correlations  at  Position  B on 
25-Degree — Rounded  Edge  . 

(60  feet  per  second,  * 8 x lO"’) 
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Figure  47b  - U^  - 100  Ft/Sec 

Figure  47  - Narrowband  Mean-Square  Pressure  Levels,  Measured  in 
62.5-Hertz  Bands  at  Various  Positions  on 
25-Degree  Knuckle  Edge 


Broadband  Space-Time  Correlations  of  Pressure  and  Velocity  above  25-Degree  Knuckle  Edge 
ssure  is  measured  at  position  F;  velocity,  above  microphone  6*  * 0.596  inch.) 


Figure  54  - Correspondence  of  Pressure-Velocity  Correlation  with 
Mean  Shear  for  Position  G on  25-Degree  Rounded  Edge 
at  60  Feet  per  Second 
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Figure  57  - Comparison  of  Autospectra  of  Pressures  on  25-Degree 
Wedge  with  Cross  Spectrum  of  Pressures  on  Opposii e Sides 

of  Trailing  Edge 

(Measurements  made  in  12.5-hertz  bands  at  El.) 


APPENDIX  A 

VARIANCE  OF  ESTIMATE  OF  CROSS  SPECTRAL  DENSITY 

In  the  fast  Fourier  transform  analyses,  the  random  variables  x(t) 
and  y(t)  are  transformed  as 


X 


i 


Xj^(t) 


2TTift 

e 


dt 


(26) 


where  X^(f)  is  the  calculated  transform  of  the  ith  sample,  and  T^^  is  the 

sample  length.  The  variables  X^(f)  and  Y^(f)  are  random  variables  and 

the  sample  of  the  cross  spectral  density  [X^(f)Y*(f)]  = [XY]^  is  also 

a random  variable;  a total  of  N Independent  samples  are  taken.  The 
ideal  cross  spectral  density  function  is 

<|)^y(f)  - < X^(f)Y*(f)  > = (XY)  (27) 


where  the  brackets  < > denote  a true  ensemble  average. 
The  estimate  of  the  cross  spectra]  density  is 

N 

“ N 2^  X^(f)Y*(f)  - (XY) 
1=1 

Now,  we  define  the  variance  of  the  sample  estimate 

0^  = < (XY  - XY)^  > 

xy 


= < 


1 N 


[(XY)^  - (XY)]‘  > 


1=1 


and  the  sample  mean  as 


(28) 


(29) 
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Expanding  Equation  (29)  we  obtain 


N 


(XY)][(XY)j  - (XY)]  > (31) 

1-1  j-1 



- I t(XY)^  - (3^)^] 

The  variables  [(XY)^  - (XY) ] and  [(XY)^  - (XY) ] are  statistically 

36 

Independent  and  so  (Cramer  (1945)  ) 

< [(XY)^  - (XY)][(XY)j  - (XY)]  > - < (XY)^  -(XY)  > < (XY)j  - (XY)  > 

- 0 

Assuming  that  the  random  variables  x and  y are  Gaussian  distributed,  the 
correlation  function  satisfies  the  relationship  (Frenklel  and  Klebanoff 

(1967)^^) 


2 2 2 2 2 2 1/2  — 2 
X Y - X Y + 2[X  Y (XY) 


(32) 


36 

Cramer,  H, , "Mathematical  Methods  of  Statistics,"  Princeton  University 
Press,  Princeton,  N.J.  (1945). 

37 

Frenklel,  F.N.  and  P.S.  Klebanoff,  "Higher  Order  Correlations  In  a 
Turbulent  Field,"  Physics  of  Fluids,  Vol.  10,  p.  507  (1967). 


Substituting  Equation  (32)  into  Equation  (31)  we  obtain 


2 1 2 2 2 2 1/2  — — 2 
^ [ X Y + 2(X  X ) ' XY  - (XY)  ] 


In  terms  of  the  spectral  density,  this  can  be  written 


fl  . 

/V'M  1 

L ^xx<^> 

Wxx<'W  J 

(33) 


We  have  assumed  that  <t>  (f)  <{)  (f).  When  x and  y are  uncorrelated, 

XX  yy 

(p  (f)  = 0 so  that  the  standard  deviation  of  the  estimate  is 
^xy 


i°<yi  ,/T 

♦„(£)  ■ V » 


(34) 


When  X and  y are  identical,  the  4>j^y(f)  = that  the  standard 

deviation  of  the  estimate  of  the  auto  spectrum  is 


a 

In  the  case  of  an  analog  cross  spectrum  analysis  N = 2AfT,  where  Af  is  the 
analysis  bandwidth,  and  T is  the  averaging  time.  Equation  (35)  is  the 
well-known  expression  for  the  error  in  the  estimate  of  an  autospectral 
density. 


»xv>  I 
xx<» 


1/^ 


(35) 
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APPENDIX  B 

DIMENSIONLESS  PROFILES  OF  MEAN  AND  TURBULENT  VELOCITIES 


Most  of  the  boundary-layer  profiles  of  mean  and  turbulent  velocity 
are  given  as  functions  of  U^  and  y/6.  The  streamline  velocity  component 
is  the  unfiltered — broadband-root-mean-squared  value.  With  only  a few 
exceptions  the  profiles  are  given  in  each  Instance  for  two  speeds 
U = 50  or  60  ft/sec  and  U = 100  feet  per  second.  Measurements  are 
given  for  all  the  beveled  trailing  edges;  see  Figures  58  through  66. 


Figure  59  - Mean  Velocity  Profiles  on  25-Degree — Rounded-Bevelcd 
Edge  at  U = 100  Feet  per  Second 
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Figure  60  - Meen  Velocity  Profiles  on  25-Degree — Rounded-Beveled 
Edge  at  U =60  Feet  per  Second 
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Figure  61  - Turbulence  Intensity  Profiles  above  25-Degree 
Rounded-Beveled  Edge  at  U =60  Feet  per  Second 
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Figure  62  - Turbulence  Intensity  Profiles  above  25-Degree- 
Rcunded-Beveled  Trailing  Edge  at  “ 100  Feet 
per  Second 
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Figure  63  - Mean  Velocity  Profiles  above  25-Degree  Knuckle 
Edge  at  = 50  Feet  per  Second 
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Figure  64  - Mean  Velocity  Profiles  above  25-Degree  Knuckle 
Edge  at  100  Feet  per  Second 
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Figure  65  - Turbulence  Intensity  Profiles  above  25-Degree — 
Knuckle-Trailing  Edge  at  “ 50  Feet  per  Second 
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Figure  66  - Turbulence  Intensity  Profiles  above  25-Degree 
Knuckle-Trailing  Edge  at  = 100  Feet  per  Second 
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APPENDIX  C 

STABILITY  ANALYSIS  OF  WAKES 

Several  characteristics  of  the  velocity  and  pressure  fields  of  the 

trailing  edges  suggest  that  the  flow  may  be  examined  using  wave 

mechanics.  The  growth  of  coherent  disturbances  in  the  wakes  are  shown 

in  Figures  16,  22,  and  23  as  well  as  the  wavelike  translation  of  velocity 

perturbations  downstream  show  the  plausibility  of  calculating  the 

stability  of  disturbances.  Similar  efforts  were  undertaken  by  Hanson^^ 

38 

for  turbulent  wakes  and  by  Crow  and  Champagne  (1971)  for  turbulent  jets. 

In  a recent  publication,  calculations  of  the  stability  of  wakes  have 

23 

been  undertaken  for  laminar  cases  by  Sato  and  Kuriki.  In  their 
experiments,  they  characterized  the  downstream  growth  of  perturbations 
in  the  streamline  velocity  in  the  laminar  near  wake  of  a flat  plate. 

The  observed  amplification  rates  and  wave  speeds  were  closely  approximated 
by  a theory  based  on  infinite  Reynolds  number.  Use  of  the  stability 

31 

theory  in  cases  Involving  turbulent  flow  has  been  attempted  by  Landahl 
who  interpreted  the  effect  of  the  turbulence  on  the  stability  as  a 
forcing  mechanism.  He  derived  a nonhomogeneous  Orr-Sommerf eld  equation 


(U(y)-c)  (4)"-k^4))  - U''4>  <({.^'^-2k^(f."+k^4.)  = q/ak 

R 

e 

for  small  disturbances  in  an  incompressible  shear  flow  of  velocity 
distribution  U(y).  The  wavelike  perturbation  velocities  have  been 
assumed  to  be  of  the  form 


Crow,  S.C.  and  F.H.  Champagne,  "Orderly  Structure  in  Jet  Turbulence," 
Journal  of  Fluid  Mechanics,  Vol.  48,  pp.  547-591  (1971). 


I 


I 


u ■ u(y)  e 


la(x-ct)+i3z 


/V  /'z  X ia(x-ct)+i0z  _ j A/  \ ioi(x-ct)-l3z 
V v(y)  e ' = - ia4)(y)  e ' ' (37) 


w = w(y)  e 


ia(x-ct)+i3z 


where  k 


is  the  wave  number  of  the  disturbance  with  wave  speed 


“ 0)/a.  Velocity  perturbations  are  assumed  to  be  much  less  than  U(y), 


The  function  q Involves  all  of  the  turbulence  terms;  the  wavelike 
dynamics  which  are  assumed  to  be  given  entirely  by  the  left  hand  s Lde  of 
Equation  (36)  are  linearly  driven  by  q.  For  the  current  discussion  all 


velocities  have  been  normalized  on  U ; all  distances,  on  a cross  stream 

o 


length  scale,*  e.g.,  b. 

The  solution  of  Equation  (36)  for  finite  Reynolds  number  is  a 
difficult  task  and  would  have  been  outside  the  scope  of  the  current  study. 


However,  Landahl  had  developed  a complete  numerical  solution  procedure 


39  40 

(Landahl  (1966)'^  and  Kaplan  (1964)  ) which  was  modified  slightly  for  this 


study  to  be  applicable  to  wake  dynamics.  The  procedure  utilizes  an 
interactive  digital  computer  program,  originally  developed  for  a time- 
sharing capability  at  Massachusetts  Institute  of  Technology.  At  the 
Center  the  program  language  was  translated  to  Fortran  as  well  as  updated 
to  handle  complex  arithmetic  more  efficiently.  Subroutines  applicable  to 
the  wake  stability  were  Included  in  the  program. 

The  computation  routine  utilizes  a Runge-Kutta  Integration  Inward 
from  a limiting  coordinate  in  the  undisturbed  free  stream,  for  Instance, 


For  discussions  In  this  appendix  all  velocities  are  assumed  to  be 
referred  to  the  local  free-stream  velocity  U . 


39, 


Landahl,  M.T.,  "A  Time-Shared  Program  System  for  the  Solution  of  the 
Stability  Problem  for  Parallel  Flows  over  Rigid  and  Flexible  Surfaces," 
Massachusetts  Institute  of  Technology  Acoustics  and  Vibrations  Laboratory 
Report  116-4  (Oct  1966). 


40., 


Kaplan,  R.E.,  "The  Stability  of  Laminar  Incompressible  Boundary  Layers 
in  the  Presence  of  Compliant  Boundaries, " Massachusetts  Institute  of 
Technology  Acoustics  and  Vibrations  Laboratory  Report  TR  116-1  (Jun  1964). 
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y = y-^  where  U(yj^)  = 1 to  the  wake  centerline,  y * 0 where  U(o)  * Ujji  ^ U(y) 
The  solution  is  controlled  by  a linear  superposition  of  the  two  functions 
(j)j^(y)  and  <j>2(y).  The  function  4>2(y)  same  as  the  viscous  solution, 

and  <|)j^(y)  is  the  inviscid,  » solution  of  the  homogeneous  Orr- 

^ 40 

Sommerfeld  equation  as  discussed  by  Kaplan.  Boundary  conditions  are 
that 

(j)j^(y)  = e ^1^  and  = e ^l\  y > 

where 

= iaN  (1-C)  + k^,  R {y}  > 0 

K 6 

e 

At  the  centerline,  it  is  stipulated  that  the  admittances  of  the  dis- 
turbances behave  as  v/p  -*•  ^ and  u/p  = 0;  p is  the  pressure  disturbance. 

That  u/p  vanishes  is  a statement  of  antisymmetry  about  y = 0.  In  the 

computation  a is  real,  and  the  wave  speed  is  complex,  C = C + iC  . 

R X 

The  mean  velocity  profile  in  the  shear  layer  U(y)  is  given  the  same 

23 

expression  as  that  used  by  Sato  and  Kurlki 

U(y)  = 1 - (1  - u^)  exp  [-P^y]  (38) 

m 

The  cross  wake  length  scale  b is  selected  so  that 


1+U 


U(b)  = 


m 


A series  of  sample  profiles  is  shown  in  Figure  67;  they  roughly  approximate 


those  that  were  mea.'.ured  at  r^  = 0.1  inch,  and  1.25,  2,  and  4 inches 


behind  the  45-deg — rounded-beveled  trailing  edge  with  * 100  feet  per 
second.  The  fifth  profile  is  given  to  indicate  dependence  on  the  factor  P. 
An  increasing  value  of  P increases  the  slope  dU(y)/dy,  and  it  produces  a 
more  narrow  wake.  Table  4 lists  a series  of  parameters  which  are 
descriptive  of  the  profiles  in  Figure  67.  The  momentum  thickness  is 


TABLE  4 - SHEAR-LAYER  LENGTH  AND  FREQUENCY  SCALES 


t 


u 

m 

P 

0 

2b 

a 

S 

cx 

m 

a 0/b 
m 

m I 

wy^ 

U 

o 

-0.1 

0.69 

0.518 

1.7 

0.52 

0.94 

0.48 

1.6 

ig 

0 

0.69 

0.62 

1.7 

0.56 

0.90 

0.56 

1.5 

mm 

0.308 

0.69 

0.752 

1.7 

0.69 

0.9 

0.68 

1.5 

1.0 

0.6 

0.69 

0. 

1.7 

0.63 

- 

- 

- 

- 

0.22 

1.20 

0.572 

1.3 

0.69 

1.1 

0.63 

1.42 

0.95 

based  on  the  full  wake  -oo  < y < oo^  and  is  the  cross  wake  distance  be- 

tween points  of  vanishing  curvature.  The  profile  for  which  P * 0.69  and 

23 

U = 0.308  also  describes  that  of  Sato  and  Kuriki. 
m 

Stability  calculations  were  made  by  checking  with  the  cases  con- 
sidered  by  Sato  and  Kuriki.  Straight  crested,  undirectional  waves  for 
which  6=0  were  assumed  for  all  calculations.  Agreement  between  their 
measured  and  current  calculated  values  of  wave  speed,  amplification  rates, 
and  disturbance  amplitude  u(y)  were  excellent.  In  the  process  of  calcu- 
lation, the  stability  diagram  of  Figure  68  was  determined.  The  magnitude 
of  the  velocity  disturbance  can  be  written 


u|  = u(y)  exp 


(39) 


thus  for  < 0 the  wake  is  stable  to  small  disturbances.  For  Reynolds 
numbers  greater  than  8.4  this  particular  profile  is  unstable;  for 


Reynolds  numbers  greater  than  300,  the  variation  of  C with  ot  appears  to 

U b ^ 

be  independent  of  N = — — - . For  the  profiles  behind  the  45-deg  rounded 

K V 

® 4 

edge,  the  Reynolds  number  based  on  b is  1.5  x 10  for  ==  = 50  feet 


per  second.  In  all  of  the  cases  calculated  no  effect  of  Reynolds  number 


was  observed  for  N > 700. 

R 

e 
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Figure  69  shows  the  calculated  intensity  eigenfunction  u(y),  compared 

with  the  measured  intensity  at  the  position  downstream  of  the  45-deg 

rounded  trailing  edge  corresponding  to  r^  = 1.25  inch.  The  absolute 

amplitudes  have  no  significance;  the  value  u(y)  expresses  the  intensity 

eigenfunction  corresponding  to  the  designated  velocity  profile.  Other 

calculations  were  performed  for  the  series  of  profiles  shown  in  Figure  67 

at  Reynolds  numbers  of  approximately  4800.  Calculated  values  of  the 

spatial  amplification  rate,  which  is  defined  in  Equation  (39),  are  shown 

in  Figure  70  for  the  profiles  in  Figure  67.  The  wave  number  a has  been 

adjusted  to  be  equivalent  to  the  dimensionless  frequency  U)y^/U^.  The 

maximum  amplification  factor  in  each  case  occurred  at  a wave  number, 

a = u)  b/C„,  where  C_  is  the  propagation  phase  velocity.  Table  4 

summarizes  the  dimensionless  wave  numbers  and  frequencies  for  each  profile. 

The  most  similar  dimensionless  wave  number  for  all  cases  is  ot  y../b  = w y,/C„. 

m-^  f m-^  f R 

The  alternative  from  cj  y^/U  is  calculated  because  it  would  seem  to  have  a 

m f o 

more  practical  significance.  The  fact  that  these  dimensionless  frequencies 

are  all  near  unity  gives  support  to  the  frequency  scaling  adopted  in 

Section  4.  The  alternative  scale  u)  9/U  can  be  seen  by  inspection  of  the 

mo 

table  to  be  not  a constant  number  over  the  range  of  profiles.  The  profile 

for  which  U = 0.6  is  shown  to  be  stable  for  all  wave  numbers.  The 
m 

correspondence  of  calculated  values  of  C„  with  the  measured  values  of 

K 

U /U  along  y = y,  in  Figure  24  is  encouraging,  and  it  raises  a question 
C ^ L 


about  the  validity  of  the  high  values  of  near  r^ 


0.8  inch. 


Finally,  we  consider  the  spatial  amplification  of  disturbances  given 


|u(y) I = exp 


a (x/b)  d 
^R 


(b) 


The  Integral  form  of  the  amplification  factor  is  made  necessary  by  the 

fact  that  C,./C_  is  a function  of  x,  i.e.,  C,,/C  (x/b).  The  form  of  the 

IK  IK 

integral  in  Equation  (40)  reflects  the  assumption  that  the  gradient 


\ 


1 


9(C^/Cj^)/9(x/b)  has  not  a strong  Influence  in  the  growth  process.  Amplifi- 
cation factors  evaluated  at  a Cj^y^/b  = 0.9  in  Figure  70  for  the  profiles 
of  Figure  67  are  shown  at  the  bottom  of  Figure  71.  For  x/b  > 6 the 
disturbances  are  lightly  damped,  and  in  reality  the  wake  should  begin  its 

decay  process.  The  evaluation  of  Equation  (40),  assuming  that  u =0.04, 

o 

is  shown  at  the  top  of  Figure  71.  Data  along  both  y = y^  and  y = y in 

Figure  23  are  shown  for  comparison.  It  can  be  seen  that  the  calculated 

growth  rate  is  approximately  half  of  that  which  was  measured.  A similar 

38 

comparison  was  drawn  in  the  case  of  Crow  and  Champagne  in  the  case  of 
jet  instability  at  high  Reynolds  number.  Perhaps  a breakdown  in  the 
linearized  theory  leading  to  the  Orr-Sommerf eld  equation  is  responsible 
for  the  underestimated  amplification  rates. 

Also,  a comparison  of  Figures  70  and  20  shows  that  although  the 
dimensionless  forms  Uiy,/U  cr  U)y,/U  (in  the  theoretical  case  we  must 

L S L 0 

accept  = U^)  align  the  peaks  of  the  disturbances,  their  bandwldths  are 
notably  dissimilar.  Thus,  physically  the  amplification  rates  are  higher, 
and  the  bandwidth  of  the  amplification  is  more  selective,  than  theory 
would  indicate. 

41 

Landahl  (1972)  has  considered  the  Influences  of  large-amplitude 
disturbances  on  the  generation  of  waves.  The  initial  disturbance  pre- 
dictable by  the  linear  theory  is  considered  to  have  reached  sufficient 
amplitude  to  initiate  a second  instability  which  in  turn  interacts  with 
the  first.  This  causes  a complex  breakdown  in  the  wake  structures  which 
could  give  rise  to  the  high-level,  selectively-amplified  disturbances 

which  are  experimentally  observed.  The  alternative  theoretical  con- 

22 

sideration  by  Abernathy  and  Kronauer  examines  the  larger  amplitude  dis- 
turbances as  being  caused  by  a nonlinear  interaction  of  the  parallel 
shear  layers  in  the  wakes.  These  interactions  culminate  in  the 
formation  of  distinct  and  localized  regions  of  vorticity.  The  vorticity 
regions  have  a spatial  organization  that  is  similar  to  that  observed  in 
vortex  streets.  A complete  theory  of  the  overall  development  of  coherent 
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wake  disturbances  would  Incorporate  both  the  linear  and  the  nonlinear 

behavior.  The  linear  calculations  as  presented  In  this  appendix  would 

be  expected  to  provide  the  growth  rate  and  wavelength  of  a preferred, 

initially  unstable  disturbance.  This  disturbance  would  be  an  initial 

value  for  the  nonlinear  problem.  The  wavelength  of  the  least  stable 

linear  disturbance  corresponds  to  the  calculated  streamline  separation  of 

22 

the  vortices  of  the  nonlinear  theory  of  Abernathy  and  Kronauer.  Their 
results  can  be  Interpreted  In  terms  of  our  wave  number  a.  That  theory  has 
shows  a decreasing  number  of  vorticlty  concentrations  with  Increasing 
values  of  a;  for  a = 0.88,  two  oppositely  phased  vorticlty  regions  are 
formed  within  a given  wavelength.  This  value  Is  very  close  to  the  values 
of  calculated  by  the  linear  theory  of  small  disturbances  as  tabulated 
In  Table  A.  The  equivalence  Indicates  a consistence  which  makes  matching 
the  linear  and  nonlinear  theories  extremely  attractive. 
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Figure  68  - Diagram  of  the  Stability  of  a Two-Dimensional 


Wake:  P =0.69  and  U =0.31 
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Figure  69  - Cross  Wake  Distribution  of  Streamline  Velocity 
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PERMANENT  TECHNICAL  VALUE.  DESIGNATED  BY  A SERIAL  REPORT  NUMBER.  | 

(2)  DEPARTMENTAL  REPORTS,  A SEMIFORMAL  SERIES,  RECORDING  INFORMA-  J 

TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR  1 

SIGNIFICANCE.  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION.  | 

(3)  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  USUALLY  INTERNAL  | 

WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS.  NUMBERED  AS  TM  SERIES  I 

REPORTS:  NOT  FOR  GENERAL  DISTRIBUTION.  £ | 


